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ABSTRACT 


This  study  addresses  the  problem  of  uncertainty  in  the 
aggregation  of  subsystem  cost  estimating  relationships  (CERs) 
for  destroyer  type  naval  ships.   The  case  of  correlations 
between  subsystems  is  discussed  with  respect  to  total  cost 
estimates,  total  cost  variance,  and  prediction  intervals. 
The  joint  generalized  least  squares  method  (JGLS)  is  devel- 
oped and  applied  to  the  ship  cost  estimating  problem.   The 
models  analyzed  were  developed  by  D.  M.  Hernon  and  R.  R. 
McCumber,  Estimation  of  Destroyer  Type  Naval  Ship  Procure- 
ment Costs,  utilizing  the  data  base  from  Resource  Management 
Corporation  Report  CR-O58,  the  Navships  Cost  Model.   One 
model  utilizes  9  subsystem  CERs  to  obtain  a  total  cost 
estimate,  while  the  second  model  utilizes  4  subsystem  CERs. 
Three  solution  methods  were  used  for  purposes  of  comparison; 
the  JGLS  method  assuming  correlated  subsystems,  the  least 
squares  method  assuming  correlated  subsystems,  and  the  least 
squares  method  assuming  independent  subsystems.   The  JGLS 
method  is  shown  to  provide  the  most  meaningful  results  for 
correlated  subsystems.   A  computer  program  and  user's  guide 
are  provided  to  conduct  JGLS  analysis. 
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I.   INTRODUCTION 

Due  to  cost  overruns,  Congressional  concern,  and  a 
continuing  need  for  better  planning  estimates,  it  is  impera- 
tive that  new  techniques  be  developed  and  old  techniques 
refined  in  order  to  obtain  better  estimates  for  major  weapon 
system  procurement  and  construction.  Along  with  these  tech- 
niques a  better  understanding  of  the  factors  and  forces  that 
determine  cost  is  required. 

The  area  of  concern  in  this  paper  is  with  Naval  ship 
procurement  and  construction  costs.   The  cost  of  a  ship  is 
a  result  of  many  factors  and  forces  all  of  which  it  would  be 
impossible  to  define  and  measure.   Thus  this  cost  is  subject 
to  variation.   When  estimating  the  cost  of  a  ship  it  is 
necessary  to  provide  both  the  estimate  itself  and  as  clear  a 
picture  as  possible  of  the  variability  surrounding  the 
estimate . 

The  general  objective  of  this  paper  is  to  analyze  the 
variability  of  ship  cost  estimates  with  respect  to  the  Joint 
Generalized  Least  Squares  (JGLS)  estimation  technique,  and 
the  Least  Squares  (LS)  regression  estimation  technique.   It 
is  hypothesized  that  the  Joint  Generalized  Least  Squares 
Estimation  technique  takes  advantage  of  additional  informa- 
tion available  and  allows  a  more  meaningful  statement  of 
the  variance  of  a  ship  cost  estimate. 


A.   NAVAL  SHIP  PROCUREMENT  AND  CONSTRUCTION  COSTS 

Naval  ships  may  be  thought  of  as  major  weapon  systems 

consisting  of  numerous  subsystems  including: 

Propulsion 

Hull 

Weapons  and  fire  control  systems 

Auxiliary 

There  are  two  basic  methods  used  to  predict  ships  costs. 
One  method  is  to  estimate  the  construction  costs  for  each 
subsystem,  which  are  then  summed  to  obtain  total  ship  cost, 
[Ref.  k   and  6].   The  other  method  is  to  estimate  total  ship 
cost  by  using  the  performance  and  operational  characteristics 
such  as,  speed  or  number  of  missile  launchers  as  explanatory 
variables,  [Ref.  2-3]. 

The  first  method  is  the  one  of  concern  in  this  paper. 
Here  the  subsystem  costs  are  estimated  by  using  the  physical 
and  performance  characteristics  of  the  ship  as  explanatory 
variables,  normally  utilizing  the  Standard  Linear  Least 
Squares  multiple  regression  technique. \    Weight  is  often  used 
as  a  primary  explanatory  variable  in  addition  to  many  other 
characteristics.   In  order  to  make  meaningful  statements 
about  the  variability  of  the  total  cost  of  the  ship  some 
aggregation  technique  must  be  employed  on  the  individual 
subsystem  relations.   The  aggregation  technique  used  depends 
upon  the  assumptions  concerning  the  correlation  of  distur- 
bances that  exist  between  the  subsystem  relations.   If  the 
assumption  of  independence  between  subsystem  relations  can 


be  made  (  a  very  restrictive  assumption)  straightforward 
linear  combinations  may  be  used  to  estimate  total  cost 
variance  and  prediction  intervals;  call  this  method  LS 
(independence).   This  would  seem  a  highly  suspect  assumption 
due  to  the  interdependencies  of  subsystems  aboard  ships. 
A  more  reasonable  assumption  would  be  that  there  are  corre- 
lations between  subsystems.   If  these  interdependencies  are 
recognized,  methods  do  exist  to  make  more  meaningful  state- 
ments about  the  uncertainty  surrounding  total  cost.   An 
extension  of  the  least  squares  method  may  be  made  by  con- 
structing an  estimated  covariance  matrix  of  the  least  squares 
estimates,  assuming  correlated  disturbances;  call  this  method 
LS  (correlated).   Finally,  if  correlations  between  subsystems 
are  assumed,  the  joint  generalized  least  squares  technique 
allows  this  information  to  be  used  in  determining  a  set  of 
joint  estimates  and  in  constructing  an  estimated  covariance 
matrix  of  the  joint  estimates.   This  assumption  is  less  re- 
strictive than  the  independence  assumption  and  in  fact 
should  provide  more  meaningful  results;  call  this  method 
JGLS.   Discussion  and  development  of  these  three  methods  is 
contained  in  Section  II,  Nature  of  the  Problem  and  Section 
III,  Methodology  of  Joint  Generalized  Least  Squares  Estimation. 


II.   NATURE  OF  THE  PROBLEM 

A.   RECENT  SHIP  COST  ESTIMATING  MODELS  AND  THE  TREATMENT  OF 
UNCERTAINTY  IN  AGGREGATE  PREDICTIONS. 

1 .   Resource  Management  Corporation  Report  CR-O58  to 
NAVSHIPS 

Resource  Management  Corporation  derived  a  statistical 

cost  model  for  NAVSHIPS  in  which  ship  construction  costs  are 

estimated  by  subsystem  and  then  aggregating  these  costs  to 

arrive  at  Basic  Contract  Cost,  and  Total  Ship  End  cost.   The 

following  subsystem  breakdown  was  used: 

1  Hull  6   Outfitting 

2  Propulsion  7  Armament 

3  Electrical  8   Design  and  Engineering 

4  Communication  and  Control   9   Construction  Services 

5  Auxiliary 

The  summation  of  these  nine  cost  categories  plus  profit  and 
overhead  was  defined  as  basic  contract  cost  and  formed  the 
nine-subsystem  cost  model.   RMC  also  developed  a  condensed 
four-subsystem  cost  model  composed  of  the  following  categories: 

1  Hull  Group         :   Hull,  Design  and  Engineering, 

Construction  Services 

2  Propulsion  Group    :   Propulsion,  Auxiliary 

3  Armament  Group     :   Armament,  Electrical, 

Communication  and  Control 

4  Outfitting  Group    :   Outfitting 

Again,  the  summation  of  these  four  cost  categories  plus 
profit  and  overhead  was  defined  as  basic  contract  cost. 


In  determining  total  end  costs  RMC  used  NAVSHIPS 
records  of  the  Shipbuilding  and  Conversion,  Navy  (SCN)  fund. 
NAVSHIPS  hardcore  cost  was  subsequently  defined  as  basic 
contract  cost  plus  miscellaneous  end  costs  and  electronics 
end  costs.   The  total  end  cost  then  became  NAVSHIPS  hardcore 
cost  plus  weapons  end  cost.   CERs  were  developed  from  the 
data  base  for  predicting  basic  contract  cost  and  miscellane- 
ous end  cost.   No  CERs  were  developed  for  electronics  or 
weapons  end  costs. 

The  models  were  developed  for  six  different  classes 
of  ships : 

1.  Aircraft  Carriers     ^ .   Auxiliary 

2.  Destroyer  5.   Amphibious 

3.  Submarine  6.   Fatrol/minesweeping 
The  area  of  interest  in  the  Hernon  and  McCumber  paper,  which 
will  be  described  next,  and  also  in  this  paper  is  with  the 
destroyer  models  developed.   Reference  6  gives  the  full 
details  of  the  RMC  study  and  results. 

Taking  now  a  point  of  departure  from  the  general  RMC 
study  development,  their  treatment  of  uncertainty  will  be 
discussed.   The  treatment  of  uncertainty  involves  basically 
two  steps.   The  first  step  is  to  develop  the  best  set  of 
CERs  from  the  data  available  based  on  a  regression  strategy 
which  outlines  the  basic  statistical  properties  to  be  used  as 
criteria  in  the  choice.   The  second  step  is  to  analyze  and 
report  as  thoroughly  as  possible  the  variance  which  surrounds 
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the  CERs  developed  and,  particularly,  their  aggregation  to 
obtain  total  cost. 

Two  points  arise  from  the  RMC  study  with  respect  to 

the  treatment  of  uncertainty.   In  the  first,  RMC  uses  the 

2 

coefficient  of  determination,  R  ,  as  a  measure  of  the  good- 
ness of  fit  of  the  regression  equation  to  the  data,  where 


*       —  ? 
£  (Y.  -Y)^ 


R2  =  ^^ 


I    (Y.  -Y)d 

i   1 


The  ratio  of  the  sum  of 
squares  explained  to  the 
total  sum  of  squares  of  Y 
adjusted  for  the  mean. 


This  statistic  tends  to  overstate  the  goodness  of  fit  as  it 
does  not  take  into  consideration  the  degrees  of  freedom.   A 
more  acceptable  statistic  would  have  been  the  adjusted 
coefficient  of  determination  (adjusted  for  degrees  of 
freedom)  where  RAnT  is  defined  below  [see  Ref.  7] 

l   e 


R2     n-k    a        n 
KADJ  '   1   z  (y  _Y)2 
n-1     a 


where  ea  =  (Ya-Ya) 


a  is  the  observation 

The  second  point  concerning  the  reporting  of  total 
cost  variance  deals  with  prediction  intervals  on  predicted 
total  cost.   The  RMC  study  implicitly  assumes  that  the  errors 
in  the  estimates  of  costs  for  the  individual  subsystems  are 
independent  of  one  another  (non-correlated),  though  they 
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fail  to  state  this  assumption  directly.   This  seems  highly 
unlikely  since  one  shipyard  produces  several  subsystems  for 
a  ship  and  the  factors  that  contribute  to  the  errors  in  the 
estimates  of  one  subsystem  may  very  well  be  the  same  or 
partially  the  same  as  those  that  contribute  to  the  errors 
in  the  estimates  of  another  subsystem.   The  allocation  of 
overhead  is  one  example.   The  factors  that  contribute  to  the 
errors  in  estimating  costs  for  the  propulsion  and  auxiliary 
subsystems  might  very  well  be  the  same  since  these  subsystems 
are  in  actuality  very  much  dependent  upon  one  another  in  the 
ship  itself. 

Though  it  is  not  reported  directly  in  the  RMC  study, 
Total  cost  variance  is  treated  as  the  summation  of  the 
variance  not  explained  in  the  individual  subsystem  CERs . 

L     2 
V  =   E   S^ 

1=1  where  V   is  total  cost 

2 
S.    is  the  variance  of 


i 


subsystem  i 

is  the  number  of 
subsystems  in  the 
model 


Prediction  intervals  on  total  cost  are  reported  in 

the  RMC  study  as  the  summation  of  the  individual  subsystem 

prediction  intervals,  where  predicted  cost  is  the  summation 
of  the  subsystem  predicted  costs. 
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where   c.  =  predicted  cost  for  subsystem  i 

Ai  =  si2(1  +  x0'(x'x)"1x0) 

the  contribution  for  the  ith  subsystem 


Acting  as  if  the  subsystems  are  independent  may  lead  to 

substantial  overstatements  of  the  confidence  one  should  have 

in  any  given  total  cost  estimate  as  this  paper  will  endeavor 

to  show. 

2.   Estimation  of  Destroyer  Type  Naval  Ship 
Procurement  Costs,  by  P.  M.  Hernon  and 
R.  R.  HcCumber  Jr. 

This  paper  is  essentially  an  extension  of  the  study 

direction  of  the  RMC  study.   The  objective  as  stated  was, 

"...  To  develop  a  model  for  the  prediction  of  total 
procurement  cost  of  destroyer  type  naval  ships  that 
increases  in  precision  as  input  data  is  refined  and, 
hopefully,  approaches  a  level  of  quality  acceptable 
in  cost-effectiveness  studies  and  eventually  for 
fiscal  planning  purposes." 

Their  approach  was  to  first  correct  numerous  defi- 
ciencies noted  in  the  RMC  study.   Because  their  concern  was 
the  prediction  of  total  cost  rather  than  the  identification 
of  basic  contract  costs  and  separate  end  cost,  data  was 
aggregated  with  the  basic  contract  cost  data  as  follows: 

a.  Electronic  end  cost  was  added  to  command  and  control 
cost 

b.  Weapons  end  cost  was  added  to  armament  cost 

c.  Miscellaneous  end  cost  was  added  to  construction 
services  cost. 
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Next  they  examined  three  models;  a  nine-subsystem  cost 
model,  a  four-subsystem  cost  model  and  finally  a  single  cost 
estimation  equation.   The  final  step  was  to  use  the  models 
developed  to  estimate  the  total  procurement  cost  of  a  de- 
stroyer type  ship  under  development  and  compare  the  predic- 
tions to  the  best  available  NAVSHIPS  estimate. 

The  study  used  as  a  primary  criteria  for  comparing 
the  prediction  values  of  these  models  the  estimate  of  total 
cost  variance  associated  with  each  model.   They  used  two 
methods  to  estimate  total  cost  variance.   The  first  method 
was  called  the  summation  method.   This  is  the  same  as  the 
treatment  used  by  RMC  in  their  study,  i.e.,  the  summation 
of  the  variance  not  explained  in  the  individual  subsystem 
CERs .   Again,  this  method  requires  the^aasumption  that  the 
errors  of  the  individual  CERs  are  independent  of  one  another. 
Hernon  and  Mc Cumber  acknowledge  the  difficulties  with  this 
assumption  and  use  this  method  to  obtain  a  minimum  total 
cost  variance  that  may  be  attained  by  utilizing  that 
particular  set  of  model  CERs. 

The  second  method  of  total  cost  variance  estimation 
involves  the  calculation  of  a  total  cost  mean  square  residual 
value  (MSR)  for  each  model.   The  following  method  was  used 
to  calculate  this  value: 


N    L  2 

MSR  =   E   [  E   (Residual  T7)] 
i=l   .1  =  1 1J 


N  -  M  -  L 


Ik 


where   N  is  the  number  of  ships 

M  is  the  number  of  variables  utilized 
in  all  CERs  of  the  model 

L  is  the  number  of  CERs  utilized  in 
the  model 


Essentially,  the  residual  values  produced  for  each  CER  for 
a  given  ship  are  summed  producing  an  aggregate  of  the  indi- 
vidual CER  residual  values.   The  total  cost  residuals  are 
squared  and  summed  for  all  observations  (ships).   If  this 
quantity  is  then  corrected  for  degrees  of  freedom,  an  esti- 
mate of  variance  is  produced  for  a  given  model.   Hernon  and 
McCumber  considered  the  MSR  method  to  produce  an  upper  bound 
on  the  total  cost  variance  estimate,  a  value  below  which 
the  estimate  of  total  cost  variance  is  expected  to  lie. 

Hernon  and  McCumber  did  not  report  prediction  inter- 
vals in  their  paper.   Their  treatment  of  total  cost  variance 
did  not  permit  utilization  of  the  standard  linear  least 
squares  methodology  for  handling  prediction  intervals.   The 
coefficient  of  determination  which  they  reported  was  the 
same  as  the  one  which  the  RMC  study  used  and  is  subject  to 
the  comments  made  earlier  concerning  it. 

3 .   RAND  Report,  "Confidence  in  Estimated  Airframe  Costs: 
Assessment  in  Aggregate  Predictions",  F.S.  Timson 
and  P.P.  Tihansky 

This  report  addresses  the  general  problem  of  confi- 
dence measures  for  multi-equation  prediction  models  utilizing 
in  the  development  specific  results  from  airframe  cost  esti- 
mation.  The  report  describes  how  prediction  intervals  can 
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be  calculated  for  sums  of  component  regressions.   The  re- 
gressions can  be  treated  as  independent  and  either  Student 
t  or  normal  statistics  used  (the  same  treatment  as  the  RMC 
study) ,  or  the  regressions  can  be  treated  as  correlated  and 
normal  statistics  used. 

A  simple  method  permits  the  determination  of  the 
degree  of  correlation  between  the  individual  regressions. 
It  is  assumed  that  there  exists  a  correlated  normal  distri- 
bution for  total  costs.   Consider  two  regressions  derived 
from  a  common  data  base  consisting  of  the  same  number  of 
observations.   Assume  the  error  terms  of  the  two  regressions 
are  distributed  bivariate  normal  with  zero  mean  then 


(eil;  e±2)   %   N(0>^} 


fl  = 


ul     u12 

2 
a21   a2 


where  i  refers  to  the  observation  number  and  Q   is  the 
variance-covariance  matrix. 


Then  if  p  is  the  correlation  between  e.,  and  e . 0  over  all 
K  ll      i2 

a  2 

observations   p  =  .   p  may  be  estimated  from  the 


sample  as 


c1  o2 


'12 


I    (Yil  -  Yil)(Yi2  -  Yi2} 

I    CYil  "  Yil»'  I    (Yi2  "  \2~> 
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When  cost  models  are  significantly  correlated  a 
method  for  treating  the  situation  is  described.   This  method 
assumes  normality  as  a  reasonably  good  approximation  to 
the  actual  probability  distribution.   The  true  variance, 
W,  of  the  sum  of  normal  costs  is  expressed  as 


W  =  Z   E   o    (1  +  X'  (X'X)  1X  ) 
i   j    J 


where  Xfl  is  a  specified  set  of  values  of 
the  independent  variables 

X  is  an  observed  n  x  K  matrix  of 
rank  K  consisting  of  values 
taken  by  K  explanatory  variables 


is  the  covariance  between 
"^   models  i  and  j  . 

An  approximation  for  which  is 


4^ 


N 

Z   e  .  e  . 


ij     N  -  K  -  1 


the  sums  of  products  of  the 
residuals  of  models  i  and  j 
over  N  observations  for 
observation  q. 


Prediction  interval  bounds  for  the  aggregate  cost  are  thus 
approximated  at  some  confidence  level  a  by 


CQ  ±  N(a)\  jl      I      s±j.(l  +  XQ'(X'X)  1   XQ) 


where  Co  the  predicted  total  cost,  is  the  sum 
of  the  individual  predicted  costs  of  the 
subsystems . 


17 


In  comparing  prediction  interval  widths  for  the  airframe 
cost  model  for  independent  normal  and  correlated  normal 
assumptions  at  a  fixed  level  of  confidence,  correlated  normal 
distributions  generally  give  the  widest  intervals. 

It  should  be  noted  that  in  this  formulation  each 
subsystem  CER  uses  exactly  the  same  explanatory  variables 
putting  a  very  limiting  constraint  on  the  use  of  this  method 
in  multi-equation  prediction  models.   In  Naval  ship  cost 
models  it  is  highly  probable  that  there  do  exist  correlations 
between  subsystem  CERs .   The  CERs  developed  in  the  RMC 
study  and  by  Hernon  and  McCumber  involve  different  explana- 
tory variables  in  the  subsystem  CERs,  hence  the  above  des- 
cribed technique  may  not  be  applied.   The  Joint  Generalized 
Least  Squares  Estimation  Technique  is  a  method  which  will 
produce  more  accurate  results  than  the  least  squares  (inde- 
pendence) technique  under  the  circumstances  of  correlated 
residuals  and  different  explanatory  variables. 

B.   THESIS  OBJECTIVES 

The  objectives  of  the  thesis  are  twofold: 

1.  To  apply  the  method  of  joint  generalized  least 
squares  to  the  ship  cost  problem. 

2.  Develop  computer  routines  to  carry  out  the  calcula- 
tions associated  with  the  joint  generalized  least  squares 
estimation  technique. 
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III.   METHODOLOGY  OF  JOINT  GENERALIZED 
LEAST  SQUARES  ESTIMATION 

Joint  generalized  least  squares  is  a  method  for  taking 
into  consideration  information  that  may  be  available  when 
one  is  interested  in  the  combination  of  several  linear  rela- 
tions which  are  related  either  because  their  coefficients 
are  partly  the  same  or  because  their  disturbances  are 
correlated.   In  the  latter  case  one  can  say  that  the  vari- 
ables that  are  neglected  in  the  various  equations  are  partly 
the  same  or  at  least  correlated.   The  statistic  used  to 
measure  the  correlation  between  relations  or  subsystem 
errors  is 


p12 


'I    (YilAl)2  l    <Yi2-^2)2 


as  was  defined  earlier  in  the  RAND  study  on  uncertainty  in 
aggregate  predictions. 

A.   MODEL  AND  ASSUMPTIONS 

A  complete  development  of  the  Joint  Generalized  Least 
Squares  Method  is  contained  in  Theil  [Ref.  7].   In  using  the 
joint  generalized  least  squares  method  it  is  desired  to 
formulate  the  estimates  of  the  parameter  vectors  of  several 
subsystem  equations  simultaneously.   Suppose  there  are  L 
equations  of  the  form 
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J     J  J     J 


j  =  1 ,  .  .  .  ,  L 


where   y.   is  a  column  vector  of  n  observations 

X.   is  an  n  x  Kj  matrix  of  values  taken 
by  the  Kj  explanatory  variables  of 
the  jth  subsystem  relation 

ft.   is  the  corresponding  parameter  vector 
J   to  be  estimated 

e .   is  a  disturbance  vector 
Suppose  also  that  each  of  these  L  subsystem  equations  meets 
the  assumptions  of  the  standard  linear  regression  model.   It 
is  desired  to  combine  the  L  equations  in  the  following  form: 


1 


X 


1 


X 


0 

h 

ci 

• 

62 

£2 

• 
• 

• 

+ 

■ 

'L_ 

h. 

_CL- 

For  notational  ease  let 


Y  = 


y- 


x  = 


X- 


x, 


£  = 


81 

el 
e2 

• 

e   = 

• 

h 

eL 

Then  the  joint  generalized  least  squares  formulation  becomes 

Y  =  X3  +  e 

1.   Assumptions  for  Joint  Generalized  Least  Squares 

a.  The  standard  linear  model  holds  for  all  of  the 
L  subsystem  relations  under  consideration. 

b.  The  disturbances  of  each  subsystem  equation  are 
homescedastic  and  uncorrelated, 
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■^jV'  ^jj1 


J    1 j  ...  5  L 


where  a.,    is  the  variance  of  the  disturbance 

th 
vector  of  the  j    equation. 

c.   Disturbances  of  different  observations  but  of 

the  same  linear  relation  are  assumed  to  be  zero, 


E(e  .e  .)  =  0 
og  nj 


a^n  j  =  1 , .  .  .  ,  L 


d.   With  respect  to  the  covariance  matrix  of  the 
disturbances  of  two  different  equations: 


E(e.s^)  = 


E(eiJe1A)   E(c1(jc2£)  . 
E(e2Jeu)   E(e2.c2£)  . 


E(£njeU)   E(£njC2£)  ' 


E(elJenl) 
E(c2Jeni) 


E(e  ,e  0) 

nj  n£ 


The  diagonal  elements  are  contemporaneous  covariances  of  the 
form  (e  .e  „ )  and  are  assumed  to  be  constant  in  the  sense  of 
being  independent  of  a.   These  will  be  denoted  by  a.«.   The 
off  diagonal  elements  correspond  to  different  observations 
and  are  assumed  to  be  zero: 


E(eaj£n*>  "  ° 


a  ?   n 


—    ±   }     ...  ,  i  ■ 


Note  we  obtain 


E(ejE') 


'n1 


j , £  -  1,  .  .  .  ,L 
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which  includes  E(e.e.)  =  a.. I  as  a  special  case. 
The  disturbances  of  the  L  linear  relations  for  the  same 
observation  are  random  drawings  from  a  multivariate  popula- 
tion with  zero  mean  and  constant  covariance  matrix.   Thus 
the  covariance  of  the  complete  vector e is: 


VAR(e) 


allI 


a21I 


"LI1 


a12I 


a22I 


aL2J 


alLJ 


^L1 


^L1 


=  Z 


where   Z  -    [a . 0 ] 


The  L  x  L  matrix  Z  is  the  covariance  matrix  of  [e 


al 


*  CaL] 


for  any  a. 

The  complete  3  vector  may  now  be  estimated  by  generalized 

least  squares.   The  Aitken  estimator  for  3  is, 


3  =  [X'(Z 


1  a  I)X]  1X'(Z"1  a  I)Y 


The  Kronecker  product  of  the  m  x  n  matrix  A  and  the 
p  x  q  matrix  B  is  defined  as  the  mp  x  nq  matrix 


A  h  B 


allB  •  •  •  alnB 


a  ,B  ...  a  B 

nl         mn 


For  further  properties  of  the  Kronecker  product  see  [Ref.  7> 
p.  305]. 
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for  which  the  covariance  matrix  is, 


VAR(B)  =  [X'  (Z  -1 


I)X] 


-1 


Since  Z  is  usually  unknown,  it  is  estimated  by  the  matrix 
of  mean  squares  and  products  of  the  least  squares  residuals 


S  = 


n 


L 


[e. 


^ 


where  e.  is  the  least 

J 
squares  residual  vector 

■f~Vi 

of  the  j    equation  and 
n  is  the  number  of 
observations. 


Replacing  Z  with  S  we  obtain  the  joint  generalized  least 
squares  estimate  of  the  parameter  vector  3 


t  c  q-1 


"J-Yt  fcj-l 


b   =  [X'  (S  "-   sa  I)X]   X'  (S 


DY 


with  the  estimator  for  the  covariance  matrix 


VAR(b.)  ='[X'(S  -1  h  I)X]  1 


Thus  we  have  the  essential  information  to  carry  out  the 
joint  generalized  least  squares  estimation  technique. 


B.   TOTAL  COST  VARIANCE 

Consider  the  estimator  for  Z: 


S  = 


n 


1 


[e1  .  .  .  eL] 
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An  estimate  for  total  cost  variance  may  be  obtained  by  taking 
the  summation  of  all  the  elements  of  S 

Var  (C)  =  i'  Si      where  i  is  a  column  vector 

of  L  ones 

This  estimate  will  be  used  in  the  analysis  of  the  ship  cost 

problem. 

C.   PREDICTION  INTERVALS 

1.   Prediction  Intervals  Utilizing  Joint 
Generalized  Least  Squares 

A  method  for  determining  a  prediction  interval 

utilizing  joint  generalized  least  squares  may  be  developed. 

It  must  be  assumed  that  the  disturbances  of  the  L  linear 

relations  for  the  same  observation  are  random  drawings  from 

a  multivariate  normal  distribution  with  zero  mean  vector  and 

L 
a  constant  covariance  matrix.   Let  C  =   E   Y.  where  C  is 

J-l   J 
total  cost;  then 

C  ~  N(  E  X*.$.,  i'Ei) 
j    J  J 

where  XK  is  a  row  of  explanatory  variables 
[X^-j  , .  .  .  ,X,jt  ]  for  the  L  subsystem  relations 

Now  if   C  =   E   Y.  where  C  and  Y.  are  predicted  costs,  then 

j-i  J  3 


C/VN(£  X^g.,  VAR(E  Xs.g.)) 


Therefore 


C-CA/N(0,  i'Ei  +  VAR(E  X-.B.)) 


2k 


If  the  estimates  utilizing  the  sample  observations  are 
i'Si  for  i'Ei  and  b.,  the  joint  generalized  least  squares 
estimate  vector,  for  3,  a  prediction  interval  may  be  derived 
as 


V- 


C  ±  ri/0/^  \/i'Si  +  X-  VAR(b,)X 


(a/2)   V  ^    •    A»  ""^MJ 


l* 


2 .   Prediction  Intervals  for  Correlated  Subsystems 
Utilizing  Least  Squares  Estimation  Techniques 

It  is  possible  to  provide  prediction  intervals  by 
single  equation  least  squares  estimation  technique  under  the 
circumstances  of  correlated  residuals  and  different  explana- 
tory variables  for  each  subsystem.   When  substantial  corre- 
lations exist  between  subsystems  in  the-  model  the  below 
described  methodology  can  be  expected  to  give  more  accurate 
statements  about  uncertainty  in  aggregate  cost  predictions 
than  would  the  normal  practice  of  disregarding  the  correla- 
tions and  without  utilizing  simultaneous  estimates. 

An  estimated  covariance  matrix  of  the  least  squares 
estimates  of  the  coefficients  is  obtained  as  will  now  be 
described.   Consider  the  L  subsystem  equations  of  the  form 

yj  =  Vj  +  ej         J-1.....1. 

Postmultiply  the  sampling  error  of  the  least  squares 

estimator,  3.,  of  3.  by  the  transpose  of  the  sampling  error 

of  the  estimator,  3-,  of  3..   Then  the  covariance  matrix 

j      j 

for  selections  i  and  j  may  be  described  as 
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EC(S1-B1)(SJ.-BJ)D  =  a   (XIX^-^X^XIX.)"1 


y  j  j 


where  c^   is  the  contemporaneous 

covariance, 

X   X 
i3  j  are  observation  matrices 

(3.-3.)  are  sampling  errors. 

This  result  will  hold  if  all  of  the  assumptions  concerning 
disturbances  that  were  made  for  the  joint  generalized  least 
squares  method  are  accepted.   These  assumptions  are  much 
less  restrictive  than  those  of  the  least  squares  (independence) 
method. 

Utilizing  this  concept  the  covariance  matrices  for 
all  combinations  of  subsystem  selections  may  be  determined. 
The  estimates  for  a. .  for  all  subsystems  are  contained  in 
S,  the  matrix  of  mean  squares  and  products  of  the  residuals. 
Letting  E.  .  =  E[  ( 3  .  -  $  .  )  ( 6  •  -  3  • )  ]  a  covariance  matrix,  denoted 
VAR($),  similar  to  the  covariance  matrix  of  the  joint  esti- 
mates may  be  constructed  in  the  form  shown  below 


Ell   E12   '  '  '   E1L 


E21   E22   *  *  *   E2L 


EL1   EL2 


E.  . 


JLL 
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If  it  is  assumed  that  the  disturbances  of  the  L 
linear  relations  for  the  same  observation  are  random  drawings 
from  a  multivariate  normal  distribution  with  zero  mean  vector 
and  a  constant  covariance  matrix,  a  method  for  determining 
a  prediction  interval  may  be  developed  paralleling  the 
prediction  interval  method  utilizing  joint  generalized  least 
squares  estimates.   The  resulting  prediction  interval 
becomes 


C  ±  n(a/2)  X/i'Si  +  X*VAR(e)x; 


where   C   is  estimate  of  total  cost 

Xx  is  a  row  vector  of  explanatory 
variables  [X#1, . . .X#L ]  for  the 
L  subsystem  relations 

VAR(3)  is  the  covariance  matrix  of  the 
least  squares  estimates. 
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IV.   DATA  BASE 

There  are  three  subsets  of  the  data  base  that  require 
discussion;  contractor  bid  cost  data,  end  cost  data,  and  the 
independent  variables.   Contractor  bid  data  was  used  for 
predicting  subsystem  costs.   This  was  a  compromise,  the 
reason  for  which  was  that  bid  cost  data  was  the  most  meaning- 
ful data  available  for  the  period  under  study  (1954-1966). 
Cost  accounting  systems  differed  greatly  among  the  various 
contractors  and  NAVSHIPS,  making  it  impossible  to  obtain 
data  on  a  uniform  level  of  aggregation  from  any  other 
sources.   Bid  prices  themselves  are  subject  to  fluctuations 
due  to  factors  in  the  shipbuilding  industry  such  as  overhead 
distribution,  workload,  level  of  expertise  at  a  particular 
shipyard.   The  source  of  the  bid  data  was  the  file  of 
NAVSHIPS  Form  4282.2,  UNIT  PRICE  ANALYSIS-BASIC  CONSTRUCTION, 
which  lists  contractor  estimates  for  the  nine  different  con- 
struction cost  groups,  subdivided  into  three  categories; 
direct  labor,  direct  material  and  overhead  costs.   The 
source  of  end  cost  data  was  the  NAVSHIPS  records  of  Ship- 
building and  Conversion,  Navy  (SCN)  fund.   There  are  forty- 
one  physical  and  performance  characteristics  for  each  ship 
which  are  contained  in  the  RMC  data  base.   These  character- 
istics are  essentially  design  parameters  such  as  maximum 
speed,  maximum  draft,  hull,  propulsion  weights,  etc.   They 
are  utilized  as  candidates  for  independent  variables  in  the 
models  developed.   Reference  4  and  reference  6  give  more 
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detailed  descriptions  of  the  entire  data  base.   Appendix  A 
gives  a  description  of  Ship's  Characteristics  used  as  ex- 
planatory variables.   Appendix  B  gives  a  description  of 
basic  contract  cost  categories. 

A.   DATA  ADJUSTMENTS 

1.  Bid  Cost  Data 

Contractor  raw  bid  data  was  adjusted  in  four  specific 
ways  to  remove  cost  variations  due  to  factors  other  than 
ship's  characteristics  as  follows: 

a.  Learning  effect  -  when  the  cost  of  ship  construc- 
tion decreases  progressively  with  each  ship  in  a  procurement 
lot. 

b.  Temporal  effect  -  which  takes  into  consideration 
variations  of  prices,  productivity  and  wage  overtime.   1965 
indices  were  used. 

c.  Installation  of  government  furnished  equipment. 

d.  Cost  of  plans  from  external  sources. 
The  order  of  adjustment  was  as  follows: 

(1)  Application  of  learning  curves. 

(2)  Adjustments  using  1965  indices. 

(3)  Addition  of  the  cost  of  GFE  and  plans. 
Details  of  these  adjustments  are  found  in  Ref.  4  and  Ref.  6. 

2.  End  Cost  Data 

End  cost  data  was  adjusted  in  much  the  same  way  that 
contractor  bid  data  was  adjusted.   Three  specific  adjustments 
were  made  to  the  raw  end  cost  data. 
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a.  Treatment  for  learning  effect. 

b.  Adjustment  for  temporal  effect  using  1965  indices 

c.  An  adjustment  for  nuclear  technology  with  propul- 
sion costs. 

Details  of  these  adjustments  are  found  in  Ref.  k   and  Ref.  6. 
The  data  base  is  published  separately  as  APPENDIX  H  of  Ref. 
4  and  is  a  CONFIDENTIAL  document. 

B.   DATA  BASE  STRATIFICATION  AND  COST  GROUP  MODELS 

Hernon  and  McCumber  noted  that  some  DLG  type  ships  had 
significantly  higher  costs  in  the  areas  of  hull,  outfitting, 
construction  services,  weapons  end  cost,  and  electronics 
end  cost.   In  addition  the  DLG  ship  was  considered  to  have 
a  different  operational  mission  than  the  smaller,  less 
expensive  destroyer  type  ships.   Thus,  two  basic  data  base 
stratification  levels  were  examined: 

1.  General  Data  Base:   DD/DDG/DE/DEG/DLG   (36  Ships) 

2.  Escort  Data  Base  :   DD/DDG/DE/DEG      (27  Ships) 
Using  both  of  the  data  bases,  CERs  were  developed  using 

two  different  methods  of  cost  disaggregation  schemes  which 
are  examined  in  this  paper.   The  coding  system  as  developed 
by  Hernon  and  McCumber  is  contained  in  Table  I  and  will  be 
used  for  easy  reference. 
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TABLE  I 
DATA  BASE  STRATIFICATION  LEVELS  AND  COST  GROUP  MODELS 


CER  CODE 


MODELS  B/C  -  9  COST  GROUP  CERS 

Hull  Cost 

Propulsion  Cost 

Electrical  Cost 

Communication  &  Control  +  Electronics 
End  Costs 

Auxiliary  Cost 

Outfitting  Cost 

Armament  +  Weapons  End  Costs 

Design  &  Engineering  Cost 

Construction  Services  +  Miscellaneous 
End  Costs 

MODELS  D/E  -  k    COST  GROUP  CER 
Base  Cost  =  Hull  +  Outfitting 


General 
Level 


D-l 


Engineering  Cost  =  Propulsion  +  Electrical  D-2 
+  Auxiliary 

Payload  Cost  =  C  &  C  +  Armament  D-3 

+  Electronics  End  Cost  +  Weapons 
End  Cost 

Construction  Cost  =  D  &  E  +  Construction   D— 4 
Services  +  Miscellaneous  End  Cost 


Escort 
Level 


B-l 

C-l 

B-2 

C-2 

B-3 

C-3 

B-4 

C-4 

B-5 

C-5 

B-6 

C-6 

B-7 

C-7 

B-8 

C-8 

B-9 

C-9 

E-l 
E-2 

E-3 

E-4 
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V.   ANALYSIS 

The  general  scheme  of  analysis  was  to  apply  the  joint 
generalized  least  squares  method  as  has  been  developed  to 
the  ship  cost  problem.   In  doing  so  computer  routines  were 
developed  to  carry  out  the  indicated  calculations  associated 
with  joint  generalized  least  squares.   Computer  routines 
were  also  developed  to  solve  the  problem  utilizing  the  least 
squares  method  in  order  to  compare  the  results  of  the  two 
methods  and  more  clearly  determine  the  nature  of  gains  by 
using  the  joint  generalized  least  squares  method. 

A.   GENERAL  COMMENTS 

1.   Computer  Programs 

A  computer  program  was  developed  in  order  to  perform 
the  joint  generalized  least  squares  analysis.   Essentially, 
a  data  set  is  read  in  and  then  transformations  to  the  data 
as  programmed  are  conducted.   The  program  then  does  the 
number  of  standard  multiple  linear  regressions  required  (one 
for  each  subsystem  CER) .   The  residuals  from  the  single 
regressions  are  saved  and  utilized  to  compute  a  matrix  of 
correlations  between  selections.   Next  the  program  does  the 
joints  generalized  least  squares  calculations;  initially  S 
is  computed  from  the  least  squares  residuals,  next  b.,  the 
parameter  estimate  vector  is  computed  along  with  the  co- 
variance  matrix  V(b.).   Finally,  the  joint  generalized 

j 

least  squares  prediction  interval  values  are  calculated. 
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Details  of  the  program  and  a  users  manual  are  contained  in 
APPENDIX  D.   The  program  was  designed  to  be  easy  to  use  for 
analysis  utilizing  the  joint  generalized  least  squares  method. 
Double  precision  was  used  to  improve  accuracy  by  reducing 
roundoff  errors.   The  program  will  handle  up  to  thirty-six 
observations  and  nine  subsystem  relations.   Each  subsystem 
relation  may  have  up  to  four  explanatory  variables  (including 
the  constant  term) .   It  should  be  noted  that  larger  dimensioned 
problems  may  be  possible  but  computer  capacity  becomes  a 
constraint  to  be  concerned  with. 

A  computer  program  was  also  developed  to  perform 
the  calculations  required  to  compute  a  least  squares  predic- 
tion interval  assuming  independence  between  subsystem  CERs . 
Essentially  a  data  set  is  read  in  and  then  transformations 
to  the  data  as  programed  are  conducted.   The  program  then 
does  the  number  of  standard  multiple  linear  regressions 
required  (one  for  each  subsystem  CER) .   Then  the  information 
required  from  each  individual  regression  is  used  in  the 
calculation  of  the  prediction  interval.   This  computer  pro- 
gram is  provided  in  Appendix  D. 

Finally,  a  computer  program  was  developed  to  perform 
the  calculations  required  to  compute  a  least  squares  predic- 
tion interval  assuming  correlated  disturbances  between  sub- 
system CERs.   Essentially,  a  data  set  is  read  in  and  trans- 
formations as  programed  performed.   Next  the  covariance 
matrix  of  the  least  squares  estimates  as  described  in 
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Section  III.C.2  is  constructed.   Finally  the  calculations 
are  carried  out  to  compute  least  squares  correlated  predic- 
tion interval  values.   This  computer  program  is  provided 
in  Appendix  D. 
2.   Criteria 

The  principle  criteria  used  for  evaluating  the  joint 
generalized  least  squares  method  were  the  reduction  in  the 
covariance  matrix  of  the  estimate,  the  prediction  intervals 
computed  and  the  mean  square  residual  values  obtained  from 
each  observation.   The  first,  a  comparison  of  the  diagonal 
elements  of  the  estimated  covariance  matrix  of  the  joint 
estimates  and  the  estimated  covariance  matrix  of  the  least 
squares  estimates  gives  an  impression  of  the  gain  obtained 
by  applying  the  joint  generalized  least  squares  method 
rather  than  the  least  squares  method.   The  second,  the  pre- 
diction intervals  computed,  gives  the  effect  of  the  assumption 
of  independent  subsystem  CER  disturbances  as  compared  to 
the  assumption  of  correlated  subsystem  CER  disturbances, 
and  also  a  feel  for  the  effect  of  using  the  joint  generalized 
least  squares  estimation  technique.   In  the  third,  mean 
squared  residual  values  (MSR)  are  obtained  from  each  obser- 
vation, as  the  difference  between  observed  and  total  cost. 
When  these  total  cost  residuals  are  squared,  summed  for  all 
observations  (ships),  and  corrected  for  the  number  of 
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observations,  a  statistic  is  obtained,  for  comparison  of  a 

2 

set  of  models. 

.  3«   Heteroscedasticity 

Hernon  and  McCumber  used  a  log  linear  model  for 

their  propulsion  and  engineering  CERs  in  their  nine  CER 

group  model  and  four  CER  group  model  respectively.   This  was 

done  in  order  to  correct  for  heteroscedasticity  encountered 

with  the  linear  model.   The  joint  generalized  least  squares 

method  will  not  allow  calculation  of  prediction  intervals 

when  individual  CERs  are  of  mixed  linear  and  log  linear  form, 

therefore,  an  alternate  method  was  used  to  correct  for 

heteroscedasticity  in  these  relationships. 

Consider  the  relation  y„  =  $n  +  3n X  -  +  e   for 

° a    0     1  al    a 

a  =  l,...,n  where  e   is  a  random  disturbance  with  zero 
expectation.   If  it  is  assumed  that  the  disturbance  variance 

is  proportional  to  the  square  of  the  independent  variable 

2       2 
then  a   =  kX  n  for  a  =  l,...,n  and  for  some  positive  constant 
a      al  '    '  ^ 


2 

Reference  is  made  to  mean  squared  residual  (MSR)  values 

in  two  contexts  in  this  paper.   MSR  values  refer  to  the  sum 
of  the  squared  total  cost  residuals  for  each  observation 
divided  by  the  number  of  observations.   MSR  values  adjusted 
for  degrees  of  freedom  refer  to  the  sum  of  the  squared  total 
cost  residuals  for  each  observation  divided  by  the  number  of 
observations  minus  the  number  of  parameters  estimated  minus 
the  number  of  CERs  utilized  in  the  model.   The  first  is  a 
statistic  for  comparison  between  models,  where  as  the  second 
is  an  upper  bound  estimate  on  total  cost  variance  as  de- 
scribed by  Hernon  and  McCumber  [Ref.  4]. 
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th 

k.   Dividing  the  values  of  the  a   observation  by  X  ,  is 

thus  equivalent  to  the  following  reformulation: 


y«  -   "°  +  6,  ♦  E« 


X  ,     X  ,     Ml    X  , 

al     al  al 


The  disturbances  e,/X, ,  .  .  .  e/X  -,  have  a  scalar  covariance 
matrix  (kl)  under  the  formulation  so  that  the  assumptions  of 
the  standard  linear  model  are  now  satisfied.   Scatter  plots 
of  residuals  versus  independent  explanatory  variables  were 
used  to  determine  the  variable  causing  the  heteroscedasticity 
in  each  relation  and  the  above  described  reformulation  was 
used  to  correct  for  it.   A  full  description  of  the  method  is 
found  in  Theil  [Ref.  7]. 

4.   Number  of  Observations 

Hernon  and  McCumber  conducted  an  analysis  of  resid- 
uals to  verify  the  normality  assumption  of  the  least  squares 
regression  technique  and  to  identify  any  significant  outliers 
This  resulted  in  the  deletion  of  outliers  in  the  development 
of  several  subsystem  CERs .   This  caused  the  observations  of 
a  particular  ship  to  be  used  to  calculate  some  subsystem 
CERs,  while  not  being  used  in  others.   Joint  generalized 
least  squares  requires  the  same  number  of  observations  for 
all  subsystem  relations.   Therefore,  all  observations  were 
used  in  the  joint  generalized  least  squares  analysis.   In 
Appendix  B  of  Ref.  h   Hernon  and  McCumber  report  the  CERs 
developed  including  outliers  in  the  data  base.   For  purposes 
of  comparison  with  joint  generalized  least  squares  results 


36 


these  CERs  were  used  instead  of  the  CERs  developed  by  the 
deletion  of  outliers.   Thus  the  general  data  base  has  thirty- 
six  observations  (ships)  and  the  escort  data  base  has  twenty- 
seven  observations. 

B.   COST  GROUP  CER  MODELS  ANALYSIS 

The  nine  cost  group  CER  model  consists  of  CERs  as  indi- 
cated and  coded  in  Table  1.   The  regression  strategy  used 
in  the  development  of  these  CERs  along  with  a  discussion  of 
the  results  is  found  in  Ref.  4.   It  was  necessary  to  use  a 
scale  factor  of  1/100  for  the  variables  ENGPAY,  PWRLD, 
ELEWGT,  C&CWGT,  PRAXWT,  LSW,  ARMWGT  and  PROLSW  (the  weight 
variables  used)  to  enable  inversion  of  the  matrices  required 
in  calculating  the  joint  generalized  least  squares  estimates. 
These  factors  are  reflected  in  the  computer  output  in  Appen- 
dix C  for  the  mine  group  CER  models,  but  have  been  removed 
in  the  results  reported  in  this  section. 

The  four  cost  group  CER  model  consists  also  of  CERs  as 
indicated  and  coded  in  Table  1.   No  scale  factor  was  re- 
quired.  In  deriving  the  condensed  model  an  attempt  was 
made  to  make  the  subsystem  categories  as  independent  as 
possible,  first  with  respect  to  the  nature  of  ship  design, 
and  secondly  with  respect  to  the  differences  reported  earlier 
in  shipyard  and  contractor  accounting  practices.   Details 
of  the  regression  strategy  and  development  of  the  subsystem 
CERs  are  contained  in  Hernon  and  McCumber, [Ref .  4]. 
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1.   Results 

The  following  tables  present  four  different  sets  of 
CERs  -  one  set  for  each  of  the  two  models  employing  the 
general  data  base,  and  one  set  for  each  of  the  same  two 
models  employing  the  reduced  escort  data  base^.   Following 
each  set  of  CERs  are  tables  containing  correlations  between 
selections  for  that  particular  model  and  selective  comparisons 
of  covariance  matrices  of  the  least  squares  and  the  joint 
generalized  least  squares  estimates.   The  complete  covariance 
matrices  are  contained  in  Appendix  C,  computer  output. 
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TABLE  IV 


BLOCK  DIAGONAL  COVARIANCE  MATRIX  COMPARISONS 
MODEL  B  -  9  CERs  -  GENERAL  DATA  BASE 


JGLS 


HULL 


.7293D-01  -.3523D-02 
.1960D-03 

PROPULSION 

.1244D  02  -.199^D  01   .1330D-03 

•  5012D  00  -.5009D-04 

.6428D-08 

ELECTRICAL 

.3025D-01  -.1964D-02  -.7572D-02 

.7613D-02  -.2417D-02 

•3206D-02 

COMMUNICATION  AND  CONTROL 

.2177D  01  -.1032D  01  -.4362D  00 

.6008D  00   .9698D-01 


.1095D  01 


AUXILIARY 

.6227D-01  -.5399D-02 
.5A25D-03 

OUTFITTING 

.1808D-01  -.4824D-03 
.1510D-04 


LS 

7392D-01  -.3578D-02 
.1991D-03 


1499D  02  -, 


2270D  01 
573*JD  00 


3377D-01  -. 


1597D-03 
9372D-02 


,238lD  01 


.1112D  01 
■6405D  00 


,1294D-03 
5755D-04 
.7885D-O8 


.9330D-02 
■3647D-02 
,42l4D-02 


,71^3D  00 
.1518D  00 
,1839D  01 


.6439D-01    -.5612D-02 
•5639D-03 


.1978D-01 


5352D-03 
.1675D-04 
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TABLE  VII 

BLOCK  DIAGONAL  COVARIANCE  MATRIX  COMPARISONS 
MODEL  C  -  9  CERs  -  ESCORT  BASE 


JGLS 


HULL. 


.M21D-01  -.2579D-02 
.17130D--03 

PROPULSION 

.5893D  01  -.55^0D  00 
.6322D-01 

ELECTRICAL 

.2340D-01  -.1152D-01  -.3137D-02 

.2160D-01  -.3461D-02 

.2137D-02 


COMPLICATION  AND  CONTROL 


.3781D  01     -.2604D  01 
.2230D  01 

AUXILIARY 

.8039D-01     -.8081D-02 
.9571D-03 

OUTFITTING 

.2786D-01  -.1039D-03 
.42O8D-06 


.5018D  00 
■.8792D  00 
.1^36D  01 


,^831D-01 


7706D  01 


LS 

-.2850D-02 
.l89^D-03 


-.7609D  00 
. 8684D-01 


.2498D-01  -.1232D-01  -.3373D-02 

.2770D-01  -.5219D-02 

.2792D-02 

.4243D  01  -.3150D  01   .1083D  01 

.3007D  01  -.1932D  01 
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TABLE 

IX 

BLOCK   DIAGONAL   COVARIANCE   MATRIX   COMPARISONS 

MODEL   D 

-    4    CERs    -    GENERAL   DATA   BASE 

JGLS 

LS 

BASE    D-l 

.1197D  00    -.3214D-04 

.1251D  00    -.3384D-04 

.1006D-07 

.1059D-07 

ENGINEERING    D-2 

• 

.1539D-06    -.1136D-04 

.9565D-05 

.2336D-06    -.1905D-04 

.3201D-04 

.17^5D-02 

-.9690D-02 

.3165D-02 

-.I898D-OI 

.1022D  00 

.1882D  00 

PAYLOAD   D-3 

.5081D  01    -.1726D-02 

•1422D  01 

.5498D  01    -.1991D-02 

.2016D  01 

.7672D-06 

-.1004D-02 

.9^95D-06 

-.l2|l|2jD-02 

.2470D  01 

.3598D  01 

CONSTRUCTION    D-4 

.7^68D  00    -.4696D-03 

.7886D-06 

.7844D  00    -.4966D-03 
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TABLE  XI 

BLOCK  DIAGONAL  COVARIANCE  MATRIX  COMPARISONS 
MODEL  E  -  4  CERs  -  ESCORT  DATA  BASE 


JGLS 


BASE   E-l 


.1056D  00  -.3761D-04 
.1456D-07 

ENGINEERING  E-2 

.1479IM)2  -.1452D-01 
.1657D  00 

PAYLOAD  E-3 

.1052D  02  -.4401D-02 
.2065D-05 

CONSTRUCTION   E-4 
.2226D  01     -.1993D-02 
.1932D-05 


.3108D  01 
•.1943D-02 
.3970D  01 

•.5692D-04 
.3906D-08 
.8972D-07 


LS 

.1108D  00    -.3962D-04 
.1533D-07 

.1533D-02    -.1513D-01 
.1727D  00 

.1179D  02    -.5108D-02 
.2472D-05 


.2374D  01     -.2120D-02 
.2050D-05 


.4255D  01 

•.2658D-02 

.5429D  01 

•.8667D-04 

.1330D-07 
.1238D-O6 
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The  following  table  indicates  correlations  of  greater 
than  ±  .4  for  the  four  models: 

TABLE  XII 


MODEL 
B 


E 


CORRELATION 

.424 

-  .405 

.402 

-  .449 

-  .418 

-  .50 

.577 

.441 
.555 

.490 


SUBSYSTEMS 

Outfitting  and  Propulsion 

Outfitting  and  C&C  + 
Electrical 

Construction  and  C&C  + 
Electrical 

Auxiliary  and  Propulsion 

Outfitting  and  Hull 

Outfitting  and  C&C  + 
Electrical 

Construction  and  C&C  + 
Electrical 

Engineering  and  Base 
Construction  and  Payload 

Construction  and  Payload 


These  correlation  estimates  do  indicate  significant  degrees 
of  correlation  between  subsystem  CERs  and  would  tend  to 
provide  evidence  against  a  hypothesis  of  independence 
between  subsystem  CERs. 

The  block  diagonal  covariance  matrix  comparisons 
show  a  reduction  in  the  variance  (diagonal)  elements  in  all 
cases  and  a  reduction  in  the  covariance  (off-diagonal)  ele- 
ments in  the  large  majority  of  cases  by  the  joint  generalized 
least  squares  estimates  as  compared  with  the  least  squares 
estimates.   This  may  be  attributed  to  the  fact  that  the  joint 


49 


generalized  least  squares  method  takes  into  consideration 
the  correlated  subsystem  disturbances  in  computing  the 
joint  estimates  whereas  the  least  squares  method  computes 
the  individual  subsystem  estimates  separately. 
2.   Total  Cost  Variance 

Total  cost  variance  for  each  model  was  computed 
utilizing  three  methods.   The  first  method  was  to  take  the 
summation  of  the  individual  subsystem  CER  standard  errors 
assuming  independent  subsystems.   The  second  method  was  to 
take  the  summation  of  the  entries  in  the  S  matrix  of  mean 
squares  and  products  of  all  sets  of  subsystem  CER  residuals, 
the  estimate  for  E.   The  third  estimate  was  obtained  by 
summing  the  residual  values  produced  by  each  CER  for  a  given 
ship  which  gives  the  difference  between  observed  and  predic- 
ted cost  for  that  ship  as  an  aggregate  of  the  individual  CER 
values.   The  total  cost  residuals  are  then  squared  and  summed 
for  all  observations  (ships)  and  divided  by  the  degrees  of 
freedom,  defined  as  the  number  of  observations  minus  the 
number  of  parameters  estimated  minus  the  number  of  subsystem 
CERs,  to  obtain  the  mean  squared  residual  value  (MSR) .   This 
method  was  used  by  Hernon  and  McCumber  [Ref .  4]  to  form  an 
upper  bound  on  total  cost  variance.   In  essence,  the  summation 
method  assuming  independence  provides  a  lower  bound  on  total 
cost  variance.   The  summation  of  the  S  matrix  elements  (S 
matrix  summation),  assuming  correlated  subsystems  provides 
a  midrange,  more  realistic  estimate,  between  the  upper  bound 
(MSR  method)  and  lower  bound  (LS  Summation  method) .   The 
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values  for  these  three  methods  are  contained  in  the  following 
table : 


TABLE 

XIII 

METHOD 

MODEL 

(LS  SUMMATION) 

(s 

MATRIX  SUMMATION) 

(MSR) 

B 

19.64 

123.9 

I63.O 

C 

16.12 

58.38 

526.0 

D 

7.86 

36.54 

58.3 

E 

7.04 

27.47 

43.6 

3«   Prediction  Interval  and  Predicted  Costs 

Three  prediction  intervals  were  obtained  for  each 
model  in  the  manner  described  in  Section  III.C,  e.g.,  least 
squares  with  independence,  least  squares  with  correlations 
and  joint  generalized  least  squares  with  correlations. 
The  estimated  total  ship  cost  was  computed  with  the  least 
squares  estimates  in  the  first  two  cases  and  the  joint  esti- 
mates in  the  third.   The  results  are  summarized  in  the 
following  tables.   There  is  a  table  for  each  model  which 
contains  the  prediction  interval  length  for  each  observation 
by  each  method.   Also  contained  in  each  table  are  the  abso- 
lute differences  between  actual  and  estimated  cost  obtained 
by  using  the  least  squares  estimates  and  the  joint  estimates. 
The  actual  prediction  intervals  are  contained  in  Appendix  C. 
The  signs  of  the  value  of  the  differences  between  actual 
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TABLE  XIV 


MODEL  B:   PREDICTION  INTERVAL  SUMMARY 


OBS 

LENGTH 

LENGTH 

LENGTH 

|A|LS 

| A| JGLS 

LS(IND) 

JGLS 

LS(CORR) 

ACTUAL- 
ESTIMATED 

ACTUAL- 
ESTIMATED 

1 

24.64 

28.76 

29.19 

10.99 

11.69 

2 

24.23 

28.54 

29.10 

3.47 

4.33 

3 

24.28 

28.55 

28.82 

3.75 

4.62 

k 

23.82 

28.37 

28.58 

4.09 

3.98 

5 

23.37 

28.12 

28.15 

1.92 

2.19 

6 

23-39 

28.12 

28.15 

1.89 

2.23 

7 

23.38 

28.11 

28.15 

.65 

•  32 

8 

23.39 

28.12 

28.15 

2.84 

3.19 

9 

23.38 

28.12 

28.15 

7.68 

7.93 

10 

23.39 

28.12 

28.15 

5.13 

5.47 

11 

23.39 

28.12 

28.15 

9.94 

9.60 

12 

23.39 

28.12 

28.15 

3.70 

4.04 

13 

23.37 

28.12 

28.15 

4.83 

5.04 

14 

23.37 

28.12 

28.14 

8.36 

8.56 

15 

23.97 

28.85 

29.02 

.46 

1.23 

16 

23.78 

28.74 

28.84 

2.80 

2.01 

17 

23.78 

28.74 

28.84 

3.37 

2.58 

18 

23.78 

28.74 

28.84 

3.78 

2.99 

19 

24.07 

29.12 

29.33 

2.15 

2.37 

20 

23.89 

28.99 

29.08 

5-33 

5.18 

21 

23.84 

28.84 

29.04 

2.31 

2.26 

22 

23.70 

28.51 

28.72 

1.49 

1.16 

23 

23.34 

28.27 

28.34 

3.43 

3.60 

24 

23.33 

28.27 

28.34 

24.15 

24.30 

25 

23.36 

28.29 

28.36 

1.08 

.78 

26 

24.22 

'28.78 

29.13 

.78 

1.53 

27 

23.79 

28.48 

28.68 

4.81 

4.22 

28 

23.56 

28.34 

28.40 

18.37 

18.39 

29 

23.56 

28.34 

28.40 

3.0 

2.99 

30 

24.56 

29.44 

29.59 

18.28 

19.39 

31 

24.12 

28.99 

29.15 

6.45 

7.29 

32 

23.85 

28.69 

28.81 

5.37 

4.76 

33 

23.85 

28.69 

28.81 

1.77 

1.16 

34 

24.48 

29.16 

29.60 

10.82 

12.26 

35 

23.19 

28.83 

28.93 

.20 

1.25 

36 

23.92 

28.84 

28.85 

2.81 

3.87 
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TABLE  XV 


MODEL  C:   PREDICTION  INTERVAL  SUMMARY 


OBS 

LENGTH 
LS(IND) 

LENGTH 
JGLS 

LENGTH 
LS(CORR) 

|A|LS 

ACTUAL- 
ESTIMATED 

|A|JGLS 
ACTUAL- 
ESTIMATED 

1 

22.02 

23.38 

23.90 

6.96 

8.0 

2 

21.73 

23.17 

23.56 

2.02 

2.79 

3 

21.76 

23.17 

23.57 

2.19 

2.93 

4 

20.71 

23.02 

23.35 

10.58 

10.05 

5 

20.92 

22.88 

22.98 

1.89 

1.62 

6 

20.92 

22.82 

22.95 

1.83 

1.59 

7 

20.91 

22.81 

22.95 

4.32 

4.10 

8 

20.92 

22.82 

22.95 

.89 

.64 

9 

20.92 

22.87 

22.98 

3.9 

4.17 

10 

20.92 

22.82 

22.95 

1.41 

1.65 

11 

20.92 

22.82 

22.95 

13.67 

13.42 

12 

20.92 

22.82 

22.95 

.03 

.22 

13 

20.91 

22.86 

22.97 

1.07 

1.31 

14 

20.91 

22.88 

22.97 

4.62 

4.86 

15 

21.30 

23.46 

23.58 

1.43 

1.18 

16 

21.29 

23.45 

23.56 

.41 

•  79 

17 

21.29 

23.45 

23.56 

.98 

1.36 

18 

21.29 

'  23.45 

23.56 

1.39 

1.77 

19 

21.38 

23.66 

23.77 

3.04 

2.35 

20 

21.38 

23.68 

23.79 

3.01 

3.48 

21 

20.87 

23.07 

23.17 

.24 

.76 

22 

21.19 

22.97 

23.35 

1.80 

1.75 

23 

21.05 

22.87 

23.17 

.60 

.95 

24 

21.04 

22.86 

23.17 

21.46 

21.77 

25 

21.30 

23.01 

23.40 

3.96 

3.58 

26 

21.28 

23.15 

23.38 

2.95 

3.31 

27 

21.11 

23.05 

23-24 

5.05 

4.44 
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TABLE  XVI 


MODEL  D:   PREDICTION  INTERVAL  SUMMARY 


OBS 

LENGTH 

LENGTH 

LENGTH 

|  A]  LS 

|A| JGLS 

LS(IND) 

JGLS 

LS(CORR) 

ACTUAL- 
ESTIMATED 

ACTUAL- 
ESTIMATED 

1 

22.49 

25.63 

25.75 

7.33 

6.70 

2 

22.37 

25.55 

25.68 

.44 

.33 

3 

22.36 

25.55 

25.68 

•  51 

.39 

il 

22.07 

25.38 

25.45 

7.01 

7.30 

5 

21.97 

25.32 

25.35 

.69 

.43 

6 

21.95 

25.32 

25.35 

.39 

•  15 

7 

21.93 

25.32 

25.36 

3.02 

2.74 

8 

21.97 

25.32 

25.35 

.54 

.81 

9 

21.99 

25.31 

25.35 

5.05 

5.34 

10 

21.97 

25.31 

25.35 

2.84 

3.10 

11 

21.97 

25.31 

25.35 

12.22 

11.97 

12 

21.97 

25.31 

25.35 

1.40 

I.67 

13 

21.99 

25.31 

25.33 

2.63 

2.86 

14 

21.97 

25.31 

25.34 

6.07 

6.31 

15 

22.07 

25.86 

25.90 

.84 

1.05 

16 

22.04 

25.85 

25.88 

1.90 

1.46 

17 

22.04 

25.85 

25.88 

2.47 

2.03 

18 

22.04 

25.85 

25.88 

2.89 

2.44 

19 

22.71 

25.82 

25.86 

2.65 

2.88 

20 

21.99 

25.80 

25.82 

4.33 

3.86 

21 

21.88 

25.58 

25.92 

.27 

.24 

22 

22.04 

25.53 

25.65 

4.81 

4.27 

23 

21.93 

25.46 

25.55 

.19 

.14 

24 

21.93 

25.46 

25.54 

20.54 

20.57 

25 

21.95 

25.46 

25.55 

4.0 

'  3.96 

26 

22.00 

.  25.43 

25.58 

1.49 

1.54 

27 

21.88 

25.42 

25.50 

4.86 

5.29 

28 

22.10 

25.41 

25.44 

19.61 

19.67 

29 

22.10 

25.41 

25.44 

1.76 

1.69 

30 

23.56 

26.27 

26.52 

18.78 

19.54 

31 

23.82 

26.23 

26.39 

9.08 

8.98 

32 

23.26 

26.16 

26.34 

5.74 

5.88 

33 

23-26 

26.16 

26.34 

2.14 

2.28 

34 

24.39 

27.04 

27.40 

4.0 

5.59 

35 

24.09 

25.87 

27.11 

5.25 

4.68 

36 

24.09 

25.87 

27.11 

2.59 

2.01 
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TABLE  XVII 


MODEL  E:   PREDICTION  INTERVAL  SUMMARY 


OBS 

LENGTH 
LS(IND) 

LENGTH 
JGLS 

LENGTH 
LS(CORR) 

|A|LS 

ACTUAL- 
ESTIMATED 

1  A  1 JGLS 
ACTUAL- 
ESTIMATED 

1 

21.39 

21.77 

22.03 

4.70 

4.59 

2 

21.15 

21.58 

21.79 

.88 

1.04 

3 

21.20 

21.55 

21.79 

.48 

.63 

4 

20.80 

21.47 

21.59 

9.35 

8.91 

5 

20.49 

21.32 

21.35 

2.07 

1.79 

6 

20.49 

21.29 

21.33 

1.68 

1.27 

7 

20.48 

21.29 

21.34 

2.74 

3.37 

8 

20.47 

21.29 

21.33 

.74 

.32 

9 

20.49 

21.32 

21.36 

3.39 

3.79 

10 

20.47 

21.29 

21.33 

1.56 

1.97 

11 

20.47 

21.29 

21.33 

13.53 

13.10 

12 

20.47 

21.29 

21.33 

.12 

.53 

13 

20.50 

21.34 

21.36 

.90 

1.30 

14 

20.50 

21.34 

21.37 

4.54 

4.93 

15 

20.88 

22.05 

22.11 

2.56 

2.29 

16 

20.88 

22.04 

22.11 

.48 

.20 

17 

20.88 

22.04 

22.11 

.07 

•  37 

18 

20.88 

22.04 

22.11 

•  50 

.78 

19 

20.59 

21.94 

22.00 

2.52 

2.16 

20 

20.59 

21.92 

21.99 

3.85 

4.24 

21 

20.13 

21.39 

21.44 

.99 

1.34 

22 

20.71 

21.59 

21.80 

3.25 

2.85 

23 

20.44 

21.39 

21.53 

.14 

.57 

24 

20.42 

21.38 

21.53 

20.77 

21.21 

25 

20.42 

21.38 

21.53 

3.65 

3-20 

26 

20.57 

21.46 

21.61 

2.72 

2.47 

27 

20.31 

21.34 

21.45 

4.71 

4.92 
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and  estimated  cost  for  both  the  least  squares  and  the  joint 
generalized  least  squares  estimates  are  the  same  in  all  cases 
and  have  been  dropped  in  the  table  for  simplification. 

The  prediction  interval  length  becomes  longer  pro- 
gressing from  least  squares  (independent)  to  joint  generalized 
least  squares  and  finally  to  least  squares  (correlated). 
Both  least  squares  (correlated)  and  joint  generalized  least 
squares  provide  intervals  which  are  significantly  longer  than 
those  produced  by  the  least  squares   (independent)  method, 
while  the  joint  generalized  least  squares  method  provides 
a  small  decrease  in  prediction  interval  length  over  the 
least  squares  (correlated)  method.   The  latter  two  methods 
represent  the  prediction  interval  length  more  accurately 
assuming  that  there  are  correlations  between  subsystem  CERs . 
Both  of  these  methods  require  approximately  the  same  degree 
of  computational  difficulty  and  thus  it  would  appear  that 
the  joint  generalized  least  squares  technique  offers  the 
best  alternative  in  accurately  stating  prediction  intervals. 

The  mean  square  residual  [MS-R]  value,  not  adjusted 
for  degrees  of  freedom,  was  computed  for  the  least  squares 
estimates  as  compared  to  actual  cost  as  was  the  same  value 
for  the  joint  estimates.   It  was  intended  to  use  these  values 
as  a  way  of  comparing  the  two  methods  to  see  which  provided 
the  best  estimates.   But  further  examination  led  to  the  con- 
clusion that  there  is  no  guarantee  in  using  this  technique 
that  one  method  will  be  consistently  better  than  the  other. 
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In  the  joint  generalized  least  squares  case  with  heterosce- 
dasticity  the  formulation  of  the  problem  is  such  that  the 
length  of  the  residual  vector,  after  a  transformation,  is 
being  minimized.   Hence  the  solution  we  obtain  is  for  the 
problem: 


Min  (Y  -  XB)(Z   @  I)(Y  -  X$) 

where  (E   x  I)  is  an  unknown  nonsingular 
matrix  which  is  estimated  by  (S~   x  I) 

and  for  which  b .  =  [X' (S_1  ©  I)X]"1X' (S_1  ©  I)Y   is  the 

joint  generalized  least  squares  estimator  as  developed  in 

Section  III.   This  estimator,  b.,  in  fact  does  not  minimize 

J 

the  length  of  the  residual  vector.   Therefore  the  MSR  method 
will  not  discriminate  between  the  two  methods  of  estimation 
in  this  case.   For  a  complete  discussion  of  this  problem  see 
Theil  [Ref.  7]. 

The  MSR  results  are  contained  in  the  following  table: 

TABLE  XVIII 
MODEL 

B       C      D      E 
LS  MSR        40. k        27-5   ^9 . 7   58.3 
JGLS  MSR      50.4    32.0    61.2    36.2 


In  models  B,  C  and  D  the  least  squares  estimates  come  out 
ahead  whereas  in  model  E  the  joint  estimates  come  out  ahead 
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The  total  cost  predictions  by  the  least  squares 
estimates  and  the  joint  estimates  are  very  close  when  com- 
pared to  each  other.   The  Summary  Tables  present  these 
values  for  comparison.   On  a  worst  case  basis  the  differences 
between  these  two  methods  for  each  model,  indicating  the 
particular  observation,  are  as  follows: 

MODEL/OBSERVATION 


B/3^      C/l      D/36      E/7 

SIZE  OP  DIFFERENCE   1.H4       1.04      1.59       .63 

$/MILLION 


In  no  case  is  this  difference  greater  than  three  percent 
of  total  estimated  cost. 

Utilizing  the  additional  information  available  due 
to  the  correlations  between  subsystems,  joint  generalized 
least  squares  estimates  should  provide  a  more  accurate 
weighting  of  the  variables  in  each  CER.   This  is  important 
in  assessing  the  contribution  to  subsystem  cost  of  a 
particular  variable.   For  example,  consider  the  CER  B-6 : 

OUTFITTING  COST  =  .436  +  .OOO386  [LSW]    (Least  Squares) 

OUTFITTING  COST  =  . 5^9  +  .000350  [LSW]    (Joint  Generalized 

Least  Squares) 

The  coefficient  of  light  ship  weight  (LSW)  reflects  the 
change  in  outfitting  cost  per  pound  for  a  change  in  1 ight 
ship  weight.   The  joint  generalized  least  squares  coefficient 
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should  provide  a  better  reflection  of  the  contribution  of 
light  ship  weight  to  outfitting  cost  for  this  CER. 
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VI .   CONCLUSIONS 

It  has  been  demonstrated  that  the  joint  generalized 
least  squares  estimation  technique,  by  taking  advantage  of 
the  additional  information  which  is  provided  when  correla- 
tions exist  between  the  residuals  of  subsystem  CERs  in  ag- 
gregate cost  estimating,  provides  more  accurate  and  meaning- 
ful statements  about  predictions,  the  total  cost  variance 
surrounding  the  predictions,  and  the  prediction  intervals. 
The  method  involves  an  increase  in  computational  difficulty 
when  compared  to  one  assuming  independence  between  subsystem 
CERs,  but  gives  a  more  accurate  statement  concerning  the 
estimates  and  their  uncertainty.   A  computer  program  has 
been  provided  which  is  readily  usable  in  conducting  joint 
generalized  least  squares  analysis  and  the  additional  time 
required  would  be  minimal. 

The  predictions  utilizing  the  joint  estimates  provide 
gains  in  a  more  accurate  weighting  of  the  variables  in  the 
subsystem  CERs.   This  allows  a  better  factor  analysis  of 
the  contributions  to  subsystem  cost  and  hence  to  total  cost. 

The  estimate  for  total  cost  variances  used  in  the  joint 
generalized  least  squares  method  provide  a  more  meaningful 
statement  taking  into  consideration  the  correlated  distur- 
bances.  The  least  squares  estimate  assuming  independence 
tends  to  substantially  understate  the  total  cost  variance 
in  the  situation  analyzed. 
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The  prediction  intervals  obtained  utilizing  the  joint 
generalized  least  squares  method  represented  a  slight  gain 
over  those  obtained  by  the  least  squares  method,  assuming 
correlated  disturbances.   Both  these  methods  provided  more 
accurate  prediction  intervals  under  the  assumptions  made \ 
however,  with  equal  computational  difficulty  the  joint 
generalized  least  squares  method  appears  the  superior  of 
the  two . 

It  should  be  emphasized  that  the  joint  generalized 
least  squares  method  involves  less  restrictive  assumptions 
than  the  least  squares  method.   The  assumption  of  correla- 
tions between  subsystem  CERs  is  a  more  reasonable  statement 
of  the  true  relationship  as  was  shown  by  the  tables  of 
correlations  between  subsystem  CERs.   The  least  squares 
method  does  in  fact  tend  to  understate  the  uncertainty 
surrounding  the  cost  estimate.   Joint  generalized  least 
squares  by  utilizing  this  additional  information  gives  more 
meaningful  results  in  the  derivation  of  the  estimate  and 
the  statement  of  uncertainty  surrounding  the  estimate. 

The  overall  gains  achieved  by  the  joint  generalized 
least  squares  method  overshadow  the  small  additional  cost. 
This  application  to  destroyer  type  ship  models  may  be  extended 
to  any  other  situation  where  one  is  interested  in  aggregate 
predictions  of  the  same  nature. 
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APPENDIX  A 


DESCRIPTION  OP  SHIP'S  CHARACTERISTICS  USED 
AS  EXPLANATORY  VARIABLES  IN  THESIS 


Characteristic   Symbol   Units 

1.   Light  Ship      LSW     long  tons 
Weight 


Definition 

Weight  of  ship  complete  with 
all  items  of  outfit,  equip- 
ment, and  machinery  but 
excluding  cargo,  stores,  crew, 
etc.  Includes  lead  ballast 
for  surface  ships  but  not 
for  submarines. 


2. 

Hull  Weight 

HULWGT 

long  tons 

3. 

Propulsion 
Weight 

PROWGT 

long  tons 

4. 

Electrical 
Weight 

ELEWGT 

long  tons 

5. 

Communication    C+CWGT 
and  Control  u/~m  nA/k 

long  tons 

6. 

Auxiliary 

Weight 

AUXWGT 

long  tons 

7. 

Outfitting 
Weight 

OUTWGT 

long  tons 

8. 

Armament 
Weight 

ARMWGT 

long  tons 

9. 

Complement 

COMP 

integer 

10. 

Maximum  Shaft 
Horsepower 

MAXSHP 

11, 


Range 


RANGE 


nautical 
miles 


These  are  the  weights  of  the 
seven  groups  as  described  in 
Bureau  of  Ships  Consolidated 
Index  of  Drawings,  Materials, 
and  Services  Related  to 
Construction  Conversion. 


Allowance  for  all  officers 
and  men  on  board. 

Total  power  that  can  be  applied 
continuously  to  the  shafts 
under  designed  operating 
conditions.  For  subs  this 
applies  to  surface  operation. 
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Characteristic    Symbol   Units 
12.   Series  SERIES    integer 


13. 

DE  or  DEG 

EEDEG 

integer 

14. 

Prototype  Dummy 

PROTO 

integer 

15. 

Number  of 
Generators 

NO-GEN 

integer 

16. 

Power  Loading 
Factor 

PWRLD 

17.   Award  Date 


18 .  Generator 
Capacity 

19.  Propulsion  and 
Auxiliary 
Weight 

20 .  Prototype 
Light  Ship 
Weight 


21.   Missile  End 
Dummy 


AWARD 


years 


TKWCPY    kilowatts 


PRAXWT    long  tons 


PROLSW 


MS-END        integer 


Definition 

A  number  that  represents 
the  position  of  a  ship  in 
a  series  of  similar  ships. 
For  example,  DD  936  to 
DDG  30  have  the  same  basic 
design  except  for  missiles 
and  constitute  a  series  of 
39  similar  ships.  DD  936 
would  be  assigned  series 
number  1  and  DDG  30  would 
be  assigned  series  number 
39. 

1  =  DE  or  DEG;  zero  otherwise 

1  =  protytype;  zero  otherwise 

Number  of  ship  service 
generators 

Ratio  of  maximum  shaft 
horsepower  to  full  dis- 
placement . 

Last  two  digits  of  the  year 
in  which  the  contract  was 
awarded . 

Total  maximum  output  of 
all  generators. 

Propulsion  weight  plus 
auxiliary  weight. 


For  prototypes,  this  takes 
on  the  value  of  LSW.  For 
non-prototypes,  it  has 
the  value  zero. 

0  =  no  launchers;  1  =  a 
launcher  at  one  end  of  the 
ship;  2  =  a  launcher  at 
each  end  of  the  ship. 
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Characteristic 

22 .  Armament 
Weight-to- 
Light  Ship 
Weight  Ratio 

23.  Engineering 
Weight 


Symbol    Units 
AR/LSW 


ENGWGT    long  tons 


2k.       Pay load  Weight    PAYWGT    long  tons 


25.   Engine  and       ENGPAY    long  tons 
Payload  Weight 


Definition 

Ratio  of  armament  weight 
to  light  ship  weight. 


Total  weight  of  all 
propulsion,  electrical, 
and  auxiliary  equipment. 

Total  weight  of  all  C+C, 
armament  and  outfitting 
equipment . 

ENGWGT  +  PAYWGT 
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APPENDIX  B 


BASIC  CONTRACT  COST  CATEGORIES 


Category 

No. 


Category  Name 
Hull  Structure 

Propulsion 


Electric  Plant 


Communication 
and  Control 


Auxiliary 


Outfit  and 

Furnishing 


Armament 


Design  and 

Engineering 

Services 

Construction 


Includes 

Shell  plating  and  planking,  logitudinal 
and  transverse  frames,  decks,  super- 
structure, armor,  etc. 

Boilers  and  energy  converters,  propul- 
sion units,  uptakes,  propulsion  control 
equipment,  feedwater  and  condensate 
system,  etc. 

Electric  power  generators,  power 
distribution  switchboards  and  cables, 
lighting  systems,  etc. 

Navigation  equipment,  interior  communi- 
cation equipment,  fire  control  systems, 
radar  systems,  radio  communications 
systems,  sonar  systems,  etc. 

Heating,  ventilating,  air  conditioning, 
plumbing,  elevators,  arresting  gears, 
rudders,  etc. 

Hull  fittings,  nonstructural  bulkheads, 
paintings,  equipment  for  work  shops, 
furnishings  for  quarters ,  etc . 

Guns  and  gun  mount,  ammunition  handling 
and  storage  systems,  other  weapon  systems 
handling  and  storage  systems,  etc. 

Contract  drawings,  working  drawings, 
technical  manuals,  lofting,  mock-up 
and  models,  etc. 

Staging,  scaffolding  and  cribbing, 
launching,  trials,  cleaning  ship, 
drydocking,  etc. 
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APPENDIX  B  (CONT.) 


END  COST  CATEGORIES 


Miscellaneous 
End  Cost 

Weapons 
End  Cost 

Electronics 
End  Cost 


Disaster  costs;  cost  of  hull,  mechanical  and 
electrical  changes;  post-delivery  costs;  etc. 

Weapons  costs  after  contractor  delivery; 
missile,  ASROC  systems;  etc. 

Electronics  costs  after  contractor  delivery; 
radar,  NTDS,  fire  control  systems;  etc. 
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APPENDIX  C 
COMPUTER  OUTPUT 
MODEL  B  -  9  CERs  -  36  OBSERVATIONS 


MLTIPLE  FEGFESSICN GEfvFAL 

SELECTICN 1 

SELECTICN  CARD  C  2  114 

LEAST  SQUARES  REGRESSION  CCEFICIENTS 


C.C4C621     C. 138546 
STC.  ERRCF  CF  ESTIMATE* 


0.6C5304 


CCVAPIANCE    MATRIX    CF    THE    ESTIMATES 

0.78264D-C1    -C.37865C-C2 
-C.37S85C-C2       C.2108CD-C3 

ADJUSTED    MULTIPLE    CCRRELATICN    COEFICIENT    0.72013C 


MULTIPLE    FEGPESSICN GENPAL 

SELECTICN 2 

SELECTICN  CARD  0  3  23331 

LEAST  SCU/ARES  REGRESSION  CCEFICIENTS 


58.66671C    -6.C91649 
STC.  EPRCP  CF  ESTIMATE* 


C.00C729 
LC13S78 


CCVARIANCE  MATRIX  CF  THE  ESTIMATES 

C.3C19CD  CC  -C.45721D-C1   0.26C75D-05 

-0.45721C-C1   C.11551D-CI  -0.  1 15S2D-05 

C.26C75C-C5  -C.11592C-C5   G.15S84D-0S 

ADJUSTED  MLTIFlc  CCRRELATICN  COEFICIENT  0.75436C 
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NULTIPLE  REGRESSION GENRAL 

SELECTION 3 

SELECTICN  CARC  C  4  21625 

LEAST  SCARES  REGRESSION  CGEFICIENTS 


C. 146642     C. 386162 
STC.  ERRCP  OF  ESTINATE= 


C. 261561 
0.29375C 


CCVARIANCE  MATRIX  OF  THE  ESTIMATES 

C.36639C-C1  -C.17416C-C3  -0.1C179D-01 
-C.174UC-C3  C.1C224D-C1  -C.29765D-02 
-0.1G179C-CL  -C.29785C-C2   0.45967D-02 

ADJUSTED  MLLTIFLE  GCRPELATICN  COEF1CIENT  0.66C927 


MULTIPLE  REGRESSION GENRAL 

SELECTION 4 

SELECTION  C/FC  C  3  21729 

LEAST  SCLARES  REGRESSION  CGEFICIENTS 


-G.9C27C1     3.173864 
STC.  ERRCP  GF  ESTIMATE= 


C. 347513 

3.643863 


CCVARIANCE  MATRIX  CF  THE  ESTIMATES 


0.259  7 2  C  CI  -C.12129C  CI 
-C.12L29C  CI  C.69671C  CC 
-0.77927C  CO   C.16562D  CC 


•0.77927D  00 
0.  165o2D  00 
C.2CC63D    01 


ACJUSTEC    NLLTIFLE    CORRELATION    CCEFICIENT    0.263167 
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MULTIPLE  PEGPESSICN GENPAL 

SELECTICN 5 

SELECTION  C£RC  C  t     118 

LEAS!  SCLARES  PEGPESSICN  CCEFICIENTS 


0.C26362     G. 224177 
STC.  ERRCF  CF  ESTIMATE= 


0.57C275 


CCVARIANCE  MATRIX  CF  THE  ESTIMATES 

0.6E175C-C1  -C.59424C-C2 
-C.59424D-C2   C.597C7D-C3 

ADJUSTED  MULTIPLE  CCRRELATICN  COEFICIENT  0.7C3E15 


f'LLTIFLE  PEGPESSICN GENPAL 

SELECTICN 6 

SELECTICN  CARD  C  7  121 

LEAST  SCUARES  REGRESSION  CGEFICIENTS 


C.43583E     C.C2E57* 
STC.  EPPCP  CF  ESTIVATE= 


0.319834 


CCVARIANCE  MATRIX  CF  THE  ESTIMATES 

G.20S47D-C1  -C.56669D-C2 

-0.56669C^C3   C.17727D-C4 

ACJLSTED  MULTIPLE  CCRPELATICN  COEFICIENT  0.7C3134 
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MLTIFLE    PEGRESSICN GENPAL 

SELEC7ICN.  .  .  ..     7 


SELECTION  C£PC  0  8  2282C 

LEAS!  SCLARES  PEGRESSICN  CCEFICIENTS 


-2.608855  9.178115 

STD.     EPRCP    CF    ESTIMATE= 


3.764265 
3.793802 


CCVARIANCE  ^TRIX  CF  ThE  ESTIMATES 


C.22551C  CI 
C.924C40  CC 


C.924C4D  CO  -0.  12052C  01 
C.271C9D  CI  -0.  13240D  01 


-C.12C52C    CI    -C.1224CC    CI       0.99298D    00 
ADJUSTED    MLTIFLE    CCRPELATICN    CCEFICIENT    0.865614 


KLLTIPLE  PEGPESSICN GENPAL 

SELECT  ICN 8 

SELECTION  C£PC  C  9  22C22 

LEAS!  StL^RES  PEGRESSICN  CCEFICIENTS 


-1.012527     C.  653961 
STC.  EPRCP  CF  ESTIMATE= 


C. 151729 
1.2C7C84 


CCVARIANCE  MATRIX  CF  THE  ESTIfATES 

C.24949C  CO  -C.885C2C-C1  -0. 127510-02 
-C.885C2D-C1  C.41244C-C1  -0.172940-03 
-0.12751C-C2  -C.172S4C-C3   C. 224200-03 

ADJLSTED  fLLTIFLE  CCRREIATICN  CCEFICIENT  0.766887 
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NILTIPLE  REGRESSION GENFAL 

SELECTICN 9 

SELECTICN    CAPC    C1C    3172223 

LEAST    SCUARES    REGRESSION    CCEFICIEMS 


3.13794C  3.949713 

STC.     ERRCP    CF    ESTIf'ATE  = 


-59.8468XS  O.C59936 

1.643151 


CCVARIANCE    MATRIX    CF    THE    ESTIMATES 

0.11594C    CI    -C.12597D    CC    -0.12864D    02  -C.4CC64D-02 

-0.12557C    CO       C.15728D    CC    -0.19555G    01  -0.12&69D-02 

-C.  12  E£  4C    C2    -C. 195550     CI       C.24C24D    03  C.52447D-C1 

-C.40C64C-02    -C.12369C-C2       0.52447D-01  C.36985D-C3 

ACJLSTED    ML7IFLE    CCRPELATICN    CCEFICIENT    0.776752 
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TABLE  CF  PREDICTION  INTERVAL  VALUES 
PI  +  PRECICTED     ACTUAL 
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TABLE    OF    LS    CORRELATED    PREDICTION    INTERVAL    VALUES 
PI    +  PREDICTED  ACTUAL  PI    - 


I 

52. 

,865238 

38 

.269414 

27, 

,270000 

23. 

,673590 

2 

46. 

,266816 

31, 

.861932 

28. 

,390000 

17, 

,457048 

3 

46, 

.5648C9 

32, 

. 156748 

28, 

,400000 

17. 

,743638 

4 

61, 

,815729 

47 

.526568 

51. 

,610000 

■3  -a 

—  —  i 

.237407 

5 

56. 

.C44C64 

41  , 

.969404 

40. 

040000 

27, 

.894744 

6 

55, 

,116949 

41 

.040287 

39, 

, 150000 

26, 

,963624 

7 

54. 

,440909 

40 

.367181 

41. 

,020000 

26, 

.293453 

8 

55. 

, 116949 

41, 

.040237 

38. 

,200000 

26, 

,963624 

9 

56, 

.332498 

42, 

.257703 

34. 

.570000 

28, 

.182909 

1C 

55, 

,116949 

41 

.040287 

35. 

,910000 

26, 

.963624 

11 

55. 

, 116949 

41, 

.040287 

5C. 

,980000 

26, 

,963624 

12 

55. 

,116949 

41 

.040287 

37. 

,340000 

26, 

,963624 

13 

56. 

,223151 

42, 

.151262 

37. 

,320000 

28, 

.079373 

14 

56. 

,126673 

42, 

.057563 

33, 

,690000 

27, 

,986254 

15 

26. 

,202169 

11  , 

.688001 

11. 

,320000 

-2, 

.326167 

16 

23, 

,  137449 

8, 

.712043 

11, 

,  510000 

-5, 

,713362 

17 

23, 

. 137449 

8, 

.712043 

12. 

,080000 

-5, 

.713362 

18 

23, 

, 137449 

8 

.712043 

12, 

,490000 

~5 

.713362 

19 

25. 

,525223 

10, 

.853086 

8. 

,700000 

-3, 

,81905  1 

20 

22, 

,184774 

7 

.605762 

12, 

,930000 

-6, 

.973249 

21 

30, 

,701544 

16 

.185230 

18. 

,490000 

1, 

.668917 

22 

40, 

,942925 

26 

.583046 

25, 

,090000 

12, 

.223167 

23 

34, 

.739514 

20, 

.566869 

23, 

,990000 

6 

.394225 

24 

34. 

-72C124 

20 

.552871 

44. 

,700000 

6, 

.385617 

25 

35, 

.518926 

21 

.338761 

20, 

,250000 

7, 

,  158596 

26 

52, 

,282515 

37 

.716865 

38. 

,490000 

23 

.151215 

27 

46, 

,427530 

32 

.083912 

27, 

,270000 

17, 

.740294 

28 

63, 

,421959 

49 

.222259 

67, 

,590000 

35, 

.022559 

29 

63, 

•43972C 

49 

.239859 

46, 

,230000 

35, 

.039998 

30 

90, 

.414635 

75 

.62C819 

93, 

,900000 

60, 

.827000 

31 

83, 

.443143 

68 

.867100 

75, 

,310000 

54 

.291057 

32 

79, 

.079522 

64 

.671006 

59. 

,300000 

50, 

.262491 

33 

79, 

.079522 

64, 

.671006 

62, 

,900000 

50, 

.262491 

34 

84, 

,022361 

69 

.223926 

80, 

,040000 

54 

.425492 

35 

73, 

,769632 

59 

.301116 

59. 

,500000 

44, 

.832599 

36 

73, 

.834053 

59 

.361273 

62, 

,  170000 

44, 

.888493 
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TABLE    CF    LS 
PI    + 


(INDEPENC)     PREDICTION    INTERVAL    VALUES 
PREDICTED  ACTUAL  PI     - 


1 
2 

"2 
4 

5 
6 
7 
8 
9 
10 
IL 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 


50 
43 
44 
59 
53 


58 

97 

29 

434250 
65 


3425 
673C 
1870 


52.73 
52.05 
.73 

.94 
.73 

.73 
52.73 
53.82 


52 
53 
52 
*.2 


53 
23 
20 
20 
20 
22 
19 
28 
38 
22 
22 
33 
49 
43 
61 
61 
87 
80 
76 
76 
81 
71 
71 


74 
67 
60 
60 
6C 
89 
55 
10 
43 
23 
21 
01 
82 
9  7 
00 
CI 
90 
92 
5  9 
59 
46 
25 
31 


5869 
0240 
5588 
0240 
4000 
C24C 
C240 
0240 
5571 
1562 
£498 
3€22 
3622 
3622 
3280 
2795 
2158 
0557 
3276 
4645 
6861 
6C96 
6547 
1038 
9516 
0284 
5568 
7549 
7549 
3455 
5260 
9580 


38. 
31. 
32. 
47. 
41  . 
41  . 
40. 
41. 
42. 
41  . 
41. 
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42. 
42. 
11  . 

8. 

8. 

8. 
10. 

7. 
16. 
26. 
20. 
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21  . 
37. 
32. 
49. 
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75. 
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59. 
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268265 
86C794 
155594 
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963307 
039195 
366074 
039185 
256592 
039185 
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039185 
15C146 
056442 
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865265 
669169 
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2220ol 
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359375 


27.2 

28.3 
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51.6 
40. 0 


39 

41. 

38. 

34. 

35. 

50. 

37. 

37. 

33. 

11. 

11. 

12. 


12.4 
8.7 
12.9 
18.4 
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23.9 
44. 
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38. 
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67. 
46. 
93. 
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59. 
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80. 
59. 
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69939 
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39993 
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39993 
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69998 
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348114 
348114 
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3C. 464706 
30.371307 
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181420 
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189982 
344112 
266452 
733521 
8  .898454 
8.389112 
9.658632 
25.605621 
20.  189301 
27.440933 
37.457642 
63.337677 
56.804947 
52.740814 
52.740814 
56.980652 
47.343135 
47.399155 


30 
29 
29 
29 


-0 

-3 

_  ^ 

-3 
-I 
-4 
4 
14 
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MODEL  C  -  9  CERs  -  27  OBSERVATIONS 


MULTIPLE  PEGFESSICK ESCCRT 

SELECTION 1 

SELECTION  C/iPC  C  2  114 

LEAS!  SCARES  REGRESSION  CCEFICIENTS 


C.5S236C     C.C97619 
STC.  EPRCR  CF  ESTIMATE= 


0.296638 


CCVARIANCE  MATRIX  CF  THE  ESTIMATES 

C.5217SC-C1  -C.3C7E9C-C2 
-C.307E9C-C2   C.2C454D-C2 

ACJLSTED  MLLTIPLE  CCRFE1ATICN  COEF-ICIENT  C. 626625 


NLLTIFLE  PECFESSICN ESCCRT 

SELECT  ICN 2 

SELECTION  CARC  0  2  132 

LEAS!  SCU4RES  REGRESSION  CCEFICIENTS 


38.C66134    -C.49215C 
STC.  ERRCF  CF  ESTIMATE= 


0.714668 


CCVARIANCE  MATRIX  CF  THE  ESTIMATES 

C.  14C26C  CC  -C.  12859C-C1 
-0.13659C-C1   C.15817C-C2 

ACJLSTED  NLLTIPLE  CORRELATION  CQEFICIENT  0.82751C 
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MJLTIPLE    PECPESSICN ESCCRT 

SELECTICN....  ..    3 

SELECTICN  CARC  C  4  21625 

LEAS!  SQUARES  REGRESSICN  COEFICIENTS 


-C.C15C67     C.6576SS 
STD.  EFRCP  CF  ESTIMATE= 


C.23C566 
0.218241 


CCVAPIANCE  M/iTR  IX  CF  THE  ESTIMATES 

C.28C98C-C1  -C.13865D-C1  -0.37943D-02 
-0.13865C-C1  C.  3  1  i  7CD-C  1  -0.587UD-02 
-0.27943C-C2  -C.58711C-C2   0  .  2  1.4C6D-02 

ACJUSTED  MLTIFLE  CCRRELATICN  CQEFICIENT  0.7390C7 


MJLTIPLE  REGRESSICN ESCCRT 

SELECTICN. 4 

SELECTICN  CAPC  0  5  21726 

LEAST  SQUARES  REGRESSION  CCEFICIENTS 


-C.445S66     3.429GC6 
STC.  ERRCR  CF  ESTIMATE= 


-1.284769 
v  3.851717 


CCVARIANCE  MATRIX  Cf    THE  ESTINATES 

£-3Z7?^£  9\    -C.35441C  CI   0.12.1890  01 

-0-3544K  CI   C.238irQ  CI  -0.21745D  01 

C.1^189C  CI  -C.21745C  CI   0.25989D  01 

ACJLSTED  MLTIFLE  CCRRELATICN  CQEFICIENT  0.C66788 
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NULTIPLE  PEC-PESSICN ESCCRT 

SELECTICK 5 

SELECTICN  C^PC  0  6  US 

LE/JS1  SCL^PES  PEGRESSICN  CCEFICIENTS 


-C.CCC657     C. 229295 
STC.  EPRCP  CF  EST1MATE= 


0.59-5721 


CCVAPIANCE  MATPIX  CF  THE  ESTIMATES 

C.1C195C  CC  -C.  1G519C-CI 
-C.1C519C-C1   C.12459C-C2 

ACJLSTED  MLTIFLE  CCRPELATICN  CGEFICIENT  0.613C91 


MLTIPLE  FEGFESSICN ESCCRT 

SELEC11CN 6 

SELECTICK  C£PC  C  7  119 

LEAST  SCL/RES  PEGRESSICN  CCEFICIENTS 


C. 571227     C.CC2215 
STC.  EPRCP  CF  ESTIMATE* 


0.253202 


CCVAPIANCE  NATPIX  CF  THE  ESTIMATES 

C.26726C-C1  -C.  14722C-C3 
-0.14722C-C3   C.5962CC-C6 

ADJLSTED  NLLTIFLt  CCRPELATICN  CCEFICIENT  0.4C143C 
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MULTIPLE  PEGPESSICN ESCCRT 

SELECT  IC^ 7 

SELECTION  C£FC  C  6  2282C 

LEAST  SCUAPES  REGRESSICN  CCEFICIENTS 


-1.CCCC13  7.991626 

STC.     EPRCP    CF    ESTIMATE^ 


2.812978 
2. 5C4982 


CCVAPIANCE  MATPIX  CF  THE  ESTIMATES 

C.12C29C  CI   C.29147C  CC  -0. 63.6120  00 

C. 291470  CC   C. 16  55  EC  CI  -0.62545D  00 

-0.63612C  CO  -C.62245D  CC   0.536160  00 

ADJUSTED  MLLTIFLE  CORRELATION  COEFICIENT  0.851478 


MULTIPLE  PEGPESSICN ESCCRT 

SELECT  ICN 8 

SELECTION  C£PC  C  9  22C22 

LEAS!  SGLAPES  PEGPESSICN  CCEFICIENTS 


-C.2C6542     C.4C6525 
STD.  ERRCF  CF  ESTIMATE^ 


0*094358 
1.2371S8 


CCVAPIANCE  MATRIX  CF  THE  ESTIMATES 

C.32421C  CO  -C. 1272*0  CC  -0.47438D-02 
-0.13724D  CC  C.7522cD-Cl  C. 965850-03 
-0.47436C-C2   C.96565D-C2   0.518580-03 

ACJLSTEC  MLLTIFLE  CCRPELATICN  COEFICIENT  0.379376 
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NLLTIPLE  FEGFESSICK ESCCRT 

SELECTICN 9 

SELECTICN  CAPC  C1C  3172223 

LEAST  SCUARES  REGRESSION  CCEFICIENTS 

3.741727     3.940985   -68.643523     0.C38777 
STD.  EPPCP  CF  ESTINATE=        1.56C89G 

CCVAPIANCE  NATPIX  CF  THE  ESTIMATES 

0.12337C    01    -C.26698D    CC    -G. 113300    02  -C.912C7D-C2 

-C.26698C    CO       C.41S1£G    GC    -0.45499D    Oi  -0.1m5590-C2 

-C.1123CC    G2    -C.454S9C     CI       0.26394D    03  C.95685D-C1 

-0.912G7OC2    -C.14559C-C2       0.95685D-G1  C.84C52D-G3 

ADJLSTED    ML7IFLE    CCRPELATICN    COEFICIENT    0.617770 
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TABLE  CF 

PREDICTICN  INI 

rERV£L  VALUES 

PI  + 

PRECICTED 

ACTUAL 

PI 

- 

1 

46.961441 

35.272428 

27.27CC00 

22, 

.583414 

2 

42.750668 

31.  168716 

28.39C000 

19, 

,586765 

2 

42.922C9C 

31  .339508 

28.400000 

19 

.756926 

4 

£2.072715 

41  .566213 

51.610000 

20, 

.05871  L 

5 

49.861508 

38.425092 

4C.C40000 

26 

.988675 

t 

48.97^£7C 

37.563356 

39. 150000 

26, 

.153842 

1 

48.322C39 

36.920174 

41.C20000 

25, 

.517308 

8 

43.97287C 

37  .563356 

38.200000 

ze 

.153842 

9 

50.  186199 

38.747780 

34.57CCCC 

27, 

.309360 

LC 

48.97287C 

37  .563356 

35.91C0C0 

26, 

,  153842 

11 

46.97267C 

37.563356 

5C. 980000 

26 

.153842 

12 

48.972670 

37.563356 

37.340000 

26, 

.  153842 

13 

5C.C639C5 

38.632108 

37.320000 

27, 

.200311 

14 

49.986179 

38.552669 

33.690000 

27, 

.117159 

If 

24.232777 

12.503649 

1 1.320000 

C, 

.774521 

16 

22. 45131C 

10.721117 

11.510000 

-1 

.CC9075 

17 

22. 45121C 

10.721117 

12. C8C0CC 

-1, 

.009075 

18 

22.451210 

10.721117 

12.490000 

-1, 

.009075 

19 

22.394977 

11  .058243 

3.7C0OO0 

-0 

.778491 

2C 

20.92  9  221 

9.094539 

12. 93C0G0 

-2, 

,750144 

21 

29.2669  17 

17.731768 

18.490000 

6, 

.  196619 

2  2 

28.334671 

26.840Co9 

25.C90000 

15, 

.345467 

22 

34.482618 

23.C49254 

23. 99CU00 

1  I, 

.615669 

24 

34.362989 

22.933478 

44.700000 

11 

.503967 

25 

25.328  19C 

23.832306 

20.250000 

12 

.326423 

26 

46.757842 

35.181267 

38.490000 

22  , 

.604691 

27 

42.244188 

31.718632 

27.270000 

20 

.193076 

MODEL  C 
JOINT  GENERALIZED  LEAST  SQUARES 
PREDICTION  INTERVAL  VALUES 


87 


gf  mrgro,rgrgrg-4.-i— * 
oooooooooo 
*  I  I  I  I  I  I 

oooooooooo 
<r  ror-  co  -J-  -*  o  <j-  cr>  o 
in  gf  r~  gf  o  r-  gjro  rg  v 
ror-rgcgcoocoaoroo^ 
r-gfrgcogogroo^-- 
-o  gr  -r  gf  — '  ro  — '— •  — *  ro 


ooo 

i  i  r 
ooo 
in*uin 
co-r 
ocrcj* 
o— io 


•tfsJO 

ooo 
I   I 

ooo 

gf-rgo 
r— rgm 
C7>00 

>t-o 


^-.  — *  rg  <t\  o  — '  —•  rg 
oooooooo 
till  I  I 
oooooooo 

o  o  cm  r-  r-  <j*  cr>  r- 

Of^comroo^-o 
so  O  o  r-  ro  m  r-  r- 
r\i  — « r-  rg  — *  o  ro  cm 


-J-ingftn 
oooo 
'III 

OOQO 
gjrgr^r- 

ino**r-<r 
rgr-gogo 
crrgp-co 
-*-*r-co 


I    I 


OQOO 

gr  <rrggo 
wo-«— o 
roin-*cr 

OCOCOvT 


mroin 
ooo 
I  f  t 
ooo 

gf  r-CO 

r-corg 
gr  coin 
inr-go 
rorgro 


ooo 
go<r  csj 

(MOJO 

rogoo 


ingj-gr 
ooo 
I  i  ' 
ooo 

rOOco 

m<Min 
rg— »in 

P--gr^ 


rgrcirorgrorgfo 
ooooooo 
I   I   I   I   I   I  t 

OOOOOOQ 

goin<j*oo-^rg 

CQOgO-«fO-^ 

rgmgorginrg-* 

<r  co  <r  — <  — *  -o  -« 


oooooooooooooooooooooooo 
I    I    II    I  I   I   I  I   I  I  I   I 

fOgf  ~*rgeg  fM  rg  — « — i  ~*  rgro  O  m  — i  — t  —(  — <  — <ro  o  O  -*  rg 

oooooooooooooooooooooooo 
I  I  I  t  I  I  I  II  I  I  I  I  I  I  I  I 
oooooooooooooooooooooooo 

mgont^cMniujogocMgjLrgjinoinmu^grr-gr^gruj 
O0Nr^mwmg3mOrOg0roC0O-^--*inrg-*grOOgrg0 

gr  in  a-  rg  o  c  -*  -*  o  cr>in  o  *■*  co  co  r-  o  ro  rg  m  cr  go  rg  go 
rg^^cgmco-*rnro— •(Mro-*gr<rgrgrrg'-*ro<rgrrg~* 


m  gr  — •  rg  oj  — *  r>  i  — •  ~^  — •  cm  ro  o  ro  ~-»  — ' -*  -*  — •  ro  o  o  — « m 
oooooooooooooooooooooooo 
I  I  I   I   I   I   I  II       I    I   I   I   I   I  I  I 

oooooooooooooooooooooooo 
ocoo-^-^cr— *~ «rO'gr-*r\jingrgrg30co-«gjcorococD 
corggrc7*m>or-roo(nm^googrr-u^r\jc>>oo^ror-r> 
gocorgcomNrrogrma)rgommgrro<rr^a3Lr,-*roin<r 
~irorgrginor-r\j— <or^r-r-gr  co— <rggOC7*g»-gr'Crrgfo 
grrgrMrocr— <r^>rro— <<rgr— 'goco^gr  gr^rocoin--'cr. 

oooooooooooooooooooooooo 
I  I   I   I   I   I   1     I   I  I   t   I 

ro.m<rmr\jrgcMrorgrg^r\jmr^rgrorgroroinrgrg»-*in 
oooooooooooooooooooooooo 

I     I     I     I     I    I     I     I     I     I     I     I     I     I     I     t     I     I     I     t     !     I     I     I 

oooooooooooooooooooooooo 

rg-rgrmr^r^r-— <iMgjrgr^a3ro-^r-v,goc>rgr».oin-j 

O  r-  o  in  rg  co  — '  r-  co  go  r—  rg  o  r*-  co  <r  -^  o  rg  *-*  o  o  ro  in 

goo^0^r-^crHigjajrgaj--mcocoingogj 

— t  i*-  in  rg  m  rg  r—  r—  — «  go  cr  rg  go  o  o*  o  o  cr*  ro  cu  a^  cr»  >o  m 

p4  r^  in  -o  ro  m  og  r-  — •  --  — '  <n  rsi  --  m  co  m  m  04  <m  m  —1  a*  in 


00 00 00000 0000 000 oooooooo 
"   1   1   1  1   I   1  t   1  r   1    1 

c\jro-^f\j(MfOfO.^-^rsif\jc^fOO'^>rrsjO"^(,r»-^^^Oro 
oooooooooooooooooooooooo 
I  I  I  I  I  I  I  I  I  I  1  I  I  I  I  I  I  I  I  I  I  I 
oooooooooooooooooooooooo 

mojuwnt- (tvu'^so^^T'^^J^^^aoi^inuio^uJirgj 
ro-jmr--*r\joao~*cum0^comcOsj'— *o^cor\ir— mr— 
mo-*<\|fMinrnr>f\jOr"i<)ro  •4--*omc\JOur>co--coco 
«j-  ro  co  tr  *-•  -4*  csi  — « ^-«  r-  — -  — •  00  o^  — « r-  -o  — •  o  ao  >r  — •  vr  -^ 

oooooooooooooooooooooooo 

I  II         I  I  I      I  I      I         I 

tMrno-^tMf^rn«-<'-*'-,rsir)roin--fn<-<(-,o>MO'-Jom 
oooooooooooooooooooooooo 

II  1  1  l  1  1  1  l  l  1  1  1  l  I  1  1  1  1 
oooooooooooooooooooooooo 
•-*  o^-^r^<-^aDtN4orMcrrnrsi'srinr-cocrr--i)-otrfO 
-^NOu^ini^^^^^^Oro^ocrr^fNjrNjfNiCTvOch-Oroin 
u>(nr*-Mp-coor,»-f\jr-u'Nrs)rsi-^o>r^inco— <-^fi-c<nrn 
^roOO>fC^^oroococoa>r--4"r--^corsjr\io-^coco 
(„nm--— *fNjaon^rg— <rsirsiNj'^-rMro-^fM— ^,c--fNjcoo,s 

oooooooooooooooooooooooo 
II    II   III    I    I   II   III 

r\j  ro  — t  rsj  -j  rsi  cm  ^1  — <  -j  og  ro  rsj  in  o  O  O  — *  r\i  »r  o  — •  -^  m 
oooooooooooooooooooooooo 
I  I  I  I  I  I  I  I  I  I  I  I  I  I  III  I  I 
oooooooooooooooooooooooo 
'T^o*-r*-m^ro^ooor-rj,ocoNTcocrp-.tMcor-r»-^-  -r 


tTsOO^'Mr—  Om— '  'TsO**4fs-ro^ro^ooor-rj,<)cONTcocrp-.r\icof^'r**^-  %; 

^^or\jf-<oorom  co-^r\jm— «r^— <o^inr,-co— -rM^^'f— «mr\icor\jri-rj<*TC 

ino-y-'rninr^-gj  Tr)-^cMgrr,tm»rin  o»rmrfi  cnin-^r  gjin^-ingrrn^c 

oo^rocuiri(7*>om  LT(>inr^NOr*--i)r\jO'£i^^%ofNj(J'<)inr--^ininin'--,-^o 


Lrc*inr—  NOr^-orMO-ci^-.ofNJC'-'v"" 

■^^.-•.^•Trn^^cNicoi^-^^mm^— **u— '— *-^rgr\j 


OOOOOOOOOOOOOOOOOOOOOOOO 

I         III      I      It         I  It      III      I 


rostrum 
OOOO 
I   1   r   1 

OOOO 

r-oco.-' 
C^coom 
OroO-« 
inm-^ 


— *  — 4  c\j  -j  r\j  rg  -j  rsj 

oooooooo 
I    I    l    l    I    I    I    l 

oooooooo 
sj-r-r— CT^o— «r\ico 
nj  o  oj>o  ro  r-  r\j  m 
ror-— <Nj-mcrcj,'M 
rgOO-^mo 
-*  rsj  in  — •  co  ro  in  o 


rM<4vr\jrsj(Mroro>r--'-^0^ 
OOOOOOOOOOOO 
I  t  I  I  1  I  I  1  I  I  I 
OOOOOOOOOOOO 

m-TvrONr^o^Tror-rsjsj- 
>tNO^-r\jr^r*-— *— <0<-^coo 
mrM-^c^rocorNjin— ^■omrNj 
ino-Of*-r-a,,incoc\j-4-r\iav 
r\j-^c^r\j-Tco-^--^--— i~*-j- 


oooooooooooooooooooooooo 
11    111   r   1  1   11   1    1 

rOst fNJrorocMro— <-*om  ^(^iino— <Oro>j-^'-j»-<-^ro 
oooooooooooooooooooooooo 
I  t  I  t  I  I  I  I  I      I  I  I  I  I   I  I   I  I      I 

oooooooooooooooooooooooo 
r\jr-cooo -or-O'-'Osj-  rsjincurj-o^comoo-r  o»r  — ' 
oooocrr\j»j-r-m  rocgu,omo----rNj--'-rr,-moco 
fNJgDro<jr\jrj*coro^^o^rsjro^-cra^p*-^-^-cor\joro— <r^ 
og^-fivomif^-o^-r-  o<Nir^sj-''Oin>j-mr-g3'-<in'Or,-o 
cOmr\ir\jror\ico--g?-  — 'cocrro  — <r\j— 'inrnr*-vrr\j.i)-^o^ 


oooooooooooooooooooooooo 


I 


I 


I 


I 


I  I 


r 


I  I 


I 


m  sj-  <N  ro  --  — *  c\j«  m  r\/  cm  o  -*  rg  >r  ~*  ro  — *  r\j  rsi  sr  — *  cm  o  ^r 
oooooooooooooooooooooooo 
I  I  I  I  I  I  I  I  1  I  I  I  I  I  I  I  I  I  I  I  I  t 
oooooooooooooooooooooooo 

fd  r-  — *  rj>  O  -T  P--  — •  *4"  -J-  >0  O  O  o*  -4-  -O  rO  r»  —•  -t  ro  O  -O  ro 
co  O  o  rg  r-  rg  (m  ro  ro  o  O  ^o  o  — •  ^O  cr  — 1  Is-  ro  o  (^  ro  Q*  r- 
O^j-  •TOxcrror-m^-'cO'-*Omgj  os<>g'4'r*'Jror"--^ro— j 
00  cj»  o  ro  >j- rg  com  <j«oni>rcoo<Mmgj^rgom— ^r^-vf 
trtn^-— *r\j— «roo^  in— •-^^-rg-^ror'i— <rg  — '— <<mct— <ni 


^-ro  •— •  rg  -*  i\trg  —•  -^  cm-*  1  grg  « 1  — '  o  o  — <— •  u^  o  — <  -^  gr 
oooooooooooooooooooooooo 
I   I   I  I    I  I   I   t   I   I  I   '   I   I  lit       I       I 

oooooooooooooooooooooooo 
c^f^OLn^gf--^C7vrMin^u^rMrocj>^*-^coro-ocoror- 
^■comLngbr^-co^Orvjin  rg(7*0*,0  0^t'r*-0>co**r-<r^jr0 
r^rgroco<^— cosj-— <0^nro--cDg3t^lnrJ,(,OLO•ocD^^J^J, 
rsj-^0Acorggj^a>cO-^r^r^r*rT,votn^3^-jrNjcr'-*-cr-- 
--p«.-jr\jni(pr.(j^uu^-  r-— j-^<"g<7*^^(MmtM— «f*-lMju'»rgin 


oooooooooooooooooooooooo 
II    III     III     I  I   I  I  I 

fM  ro  O  -*  ro  ro  m  oj  rvj  ro  — •  r\j  -^  gr  rg  ro  rg  —<  rg  m  sj-  in  ro  co 
oooooooooooooooooooooooo 
II  I  I  I  I  I  I  I  I  t  I  I  I  I  I  t  I  I  I  t  1  I 
oooooooooooooooooooooooo 
r-r-gr  — ic?>rgm— « cr- oo  c  co  co  sy  inOr--^r-r-r'--jrggo 
-^royg-cNj-^iAO'-JNfO^-minumiAcog'^g-in'-* 
(M<fOcoo^inrj'aDLnrMCorMn~icocogDrg-^r^g3crxoinro 
f*»^g3Noro—«cgrgrgcoo*iri--m  coor— gocgcoinu^-org 
•-<— ^r-ao--«-^gjrorg>rrgro— «g*  rgrj— -^crgocggDrggD 


oooooooooooooooooooooooo 
I    f  I     III    I    I   I  I   I  I  I 

gfinrOst^^r^fO(^^-gfinro<>inu^ining3r-ro>Tror«-  insor- 

000000000000000000000000  000 

1   1    1   1    t    1    t    1    1    1    t    1    t    1    1    1    1    1    1   1    1    r   1    t  111 

oooooooooooooooooooooooo  000 

rg-jfO-oCT^mgr^gr-O'M— 'rorsicoco<N<-«r-ccrrn-0'J"  OO-J" 

vrsO-Oro— 'g^r-omr-gDcoroou^in— 1  cr> go  rg t^  ro  o  co  r-00 

mm— <aDor\iLninoi^cogr%orsjcoaNCTv---'CD-ocrro--*-J  gj«j-g^ 

gDr^Oin0  0ogfcosrmrgromror:(c7'r-<-oor-cug3  _i„f^ 

mro«J"-4v-J'— '— *0<t  —*rOs3-— 'm^--^-  — «o— *-^g3r\i'co  cm— im 


co-t  grmrorgrMingo^r*g?-rorororororMrg-^ro 
000000000000000000000 
I  I  l  I  l  1  I  I  1  I  I  I  1  1  l  I  l  I  I  I  I 
OOOOOOOOOOOO OOQOOQOCO 
gD  ro  vj- -- CO  O  r*  0  CO  — *  gr  r- -^  g^  n  1  r- ^r  m  rg  0^  O 
— tr-o^c^coocor^cocorocoOinorgcr'OOO 
ro  --  rg  o  ^r  go  *v  cr-  rg  cm  --  cr  r-  o  ro  r-  — <  go  gr  in  -o 
rggro^rorogOr-g?,ro--gj— *Ocococo^-r-(M--^ 
g3fogrmcj'--rMsfmr\jcomc7>(MO,-^-*r'-— <ajr- 


OOOOOOOOOOOOO00000000000 

11    1  1   1  1  i   1   1  1    1    r  1 


oooooooooooooooooooooooo 
I  I   I   I   I    III    I     III 


-jcg-*oogrg  j*— *— * 
000000000 
II     Mill 

OOOOOOOOO 

g0  eg  r\j  gf-  — *  r-  g^  O  rg 
gOfNJgOOOcOgO^Jro 
gj- O  O  O  <f  O  — *  rg  r- 
cgOr-  -0  O  O  mrg  0s 
^^p-r-— '— <mrg— < 


rgo— <^rO— *rg^^o— J-g,ro<j-rgrfc- 
OOOOOOOOOOOOOOO 

1   I   1   t   1  1       i   r  1   1   1   1 

oooooooo 00 00 ooo 
r-co^p^fnoco-"a>g,N-jco-og" 
h-inrMingoi^orginor— cg-*goo 
rgo^rg-^roro-^r-  fOgOcointMgo 
rg—'  om  ocM-oin  o— «r-  crr---o  gr 

gf  rgrocu^  — irg-^—«  cor-  -^m  rgm 


— iro— '-^(Mrsim  o  ofMrg 
OOOO 0000 000 

I    I     I    I     I    I    I  t    I 

OOO OOOO OOOO 

^cor*-coOincoinmingo 
o  >r  m  — .  0  -4-  o  go  go  r-  ro 
^"r-rgromro— 'incoy^rg 
mrgm— *ooinor~-— ^g0<7^ 
~*V  oo^^rgrgrg--— -rorj* 


I    I 


1    I 


Or-—min 
i^-inrocr  0 

sJ-COgO 

Omror-<r 
— <ro-^rgro 


rgrorom— < 
OOOOO 
I  I  I  I  I 
OOOOO 

oro  <r— <  in 
rgorggr  go 
rOCMgrrggg 

rgCT  rc\^Q  lT 

-^-rtMrgrg 


_-_.g> 

OOO 

I  t  I 

OOO 
goro— ' 
rg^rco 
goinog 
g3  0>-. 
—i  gO  CM 


oooooooooooooooooooooooo 

I  I  I  I     I  I  I  I     I     I  I  II 

O  rg  rg  ro  o  O  — *  — «  -j  — <  o  — *  — *  ro  — '  — >  — '  O  O  <N  — •  — '  ro  — « 

OOOOOOOOOOOOOOOOOOOOOOOO 

I    I    t  I  lit  I  I 

oooooooooooooooooooooooo 
rgmgO(Ngjrgma30<0^<^cci^ro<^4*<r^o-romo,,gfON 
r-crgomu^coror-crroogrgfgO'-<o-g*roLr---jmcoo 
gr-g,-rginroinr-mgr*-^ocom-^cM— <<rinmromr— <rm 

ro  r>  ro  go  r-  rg  o  rg  rg  r-  o  — *  V*  co  O  r-  —*  CO  CO  CM  go  CO  <M  — ' 
-4  co  rg  cm  -*  -*  (S-  --  rg  ro  — *  —  O  <M  fM  — *  ^  co  gr  gO  0^  ro  cm  co 


oooooooooooooooooooooooo 

III  I  II  I     I  I    I  I 

rgrorMrogr<ru^gj-gruirgroromrogrrororomrgrorggo 
o  o  o  o  o  o  o  o  o  o  t_»  o  o  o  o  0  o  o  o  o  o  o  o  o 
I  1  l  t  I  I  1  1  l  1  1  t  l  1  1  l  1  1  l  1  l  1  1  I 
oooooooooooooooooooooooo 
(^cu'gr-r^gr^rcof^nirj^nico— -r^r-gjgjajiNg^— «ino 
ororgmogor— oggo  grgo  g*-grgocoo— ^gjrgrjMMr—u^o 
inc/*0"T  -r  coococrr-coogr*-  mcMgor,~iro-^r«-cr'Co>-gvgr 
go  coo— <co>c^r^romgrtMmrgr---<rc>roor\irgrg(7*gr 
rg  — <  — *  — '  m  ro  r-  cm  — <  gr  — >  — ^  rj*  ro  r—  m  ro  in  r  1  — » -^  m  cu  — * 


oooooooooooooooooooooooo 
I    I   t   t   I  I    I   I  I   t  I  I 

rgrO-^rg-^rgcMnirorocMcMtMinrggf  rororo^JCMro-^gO 
oooooooooooooooooooooooo 
I  I  I  I  I  I  1  t  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 
oooooooooooooooooooooooo 

gjniViajocor-i'gr--gjin'<-'0'MgriNJc>u>c\ig-'o<Mcocu 
gr^r  cmo^  gjfO'Ni-—  Lrgra^rnu'oo^tr-^gogrgjcogo^r- 
(O^fMo^rMOrgcooro^roLngrcfroromrMorgco^incocM 
rgm  Ouogr-Or\jr-o~*r-— «orgr^r*gocogorogrgr-^ro 
rg-^nro  ^-gocMgr^o-^fJ,--— *g^-^C^C^,tM--^r-^m-^u^ 


oooooooooooooooooooooooo 
I   t  I  I    II    till    II 


rgmgr  u 
oooc 
I   I   I 

oooc 


j*r\jr-  u 
r-mu"\»- 


g  -<  rg  o  O  — '  rg  ro  -*  gr  — *  — •  — •  — '  — <  rO  O  o  --  rg 

^0000000000000000000 

t  1     I  I  1  I  I  1  I  I  I  I  t       I 

^ 00 00 00 000 oooo 00 00 00 

_jr— orMn*o"»inrggj  mcu  (o^*^l^^Ja)f,l--^y,,',' 
lL^oro— •ogogjrgmgrinu^gjco'-'-g'fMmo 
jgomroooomgOrocDfOrogrr--- 'r-^f-g 
■*grrjvO''Omgogr-or-— *>u— '— '--inrgmcDrg 
r-~- -^^mgrgogrm-^goin  go— ^rg*— •— «cm  ror^-- 


oooooooooooooooooooooooo 
III    I    I    I   I  I  I      II 

-j  rg  — *  rg  rO  ro  ro  ro  ro  gf  — ■  rg  — •  -T  — « rO  rg  rg  rg  <j-  — *  rg  O  m 

oooooooooooooooooooooooo 
1  1  t  1  1  t  1  1  1  t  t  1  1  1  1  1  r  t  1  i  1  1   1 

oooooooooooooooooooooooo 
Ocr»gjr-ror-cMoingrgrgogrrgu,'*vj<r--r,i^-'Ofj'rggj 
muo- *oor-ocoromrogrO  gr  <r  ocu  — >mm  error-  r- 
niu^grrNjou^gogOOrooi^grinr-rg— ti^rnrgmr^-grgo 
cu  co  rgr- co  0-4— jgrf^-rgogm  gorgrg(jig>incr-^trro^^ 
■grcg— '— term  -ggr  <m  gjegrg— •  in— *  corny  <r  -<rjN- «cm 


00000(_>OOOOOOOOOOt->CJOOOOOO 

III  III  I       I    t       l    I    l  1 

-^rgo— 10— <^rgrgrg-^rgrggr  — trOrgrgrgin— <rgom 

oooooooooooooooooooooooo 

II    t    I  I  I  I  I  I  t    I  I  I  t  I  I  I  I  I  I    I 

oooooooooooooooooooooooo 

■grtrcoO  go<Mrg  — r«)u-\r-u^gogom<rr-rggrcugomroo 
rogoinoorgr-in  gjco(o>'a*rogom(Mcoro-g>rMmg-»oco 

OCOOCO-*U,,rMrggjomrorgCDmcggru\imin^OOO,v 
rsjrg— 'crrninc^r-mo^r-o^mr-rg-^corogo^gjogr 
rg— <rgrg—  m-—  gj-rg— 'OQ^^ro-^CDCorg-Mfrj^-gr^gr 


'w'O'— 'O  tot-*  OOO  00  00  <— 'OOOOOOOOOO 

I   I  I   I  I   I  I    I  1    I    I   *  I 

"-•rgrO-T— •-^rgoOO-^fM--r0O-JOO— *rg--0--rg 

oooooooooooooooooooooooo 
I  I  I  1  I  I  I      I  f  I  I   I    t  I 

oooooooooooooooooooooooo 
o^^r^ror^r^coror-goomcogo-rr-gjgr^o^i^mm 
y»  rvj  r\j  •— 1  rj^  gj  o  cr*  ^r  co  m  cu  go  cr  r-  r—  r—  cr  ro>  o  o  gr  ^  cr 
u^CDtrr-^f1)- *OfO— ii_rgjuxgjO,g"rr,rggomrwii>gj 
^rgo^inmrMingrco^rogr-g>cj,'Oi7ymmococMOrggor- 
iM---^rgcg-^rocosr  ^----*rg-*cMrg— *--gr  — *— 'rgcr>r- 


OOOOOOOOOOOOOOO000000000 
I   I   I  I   I  l    I  l    I       lit     1 


oooooooooooooooooooooooo 
I   I   I  I  I     I  I   I    I     II 


oooooooooooooooooooooooo 
I  I    II   II   II     I  I  '  t  '  » 


88 


TAELE    OF    LS    CORRELATED    PREDICTION    INTERVAL    VALUES 
PI    +  PREDICTED  ACTUAL  PI    - 


I 

46, 

,179220 

34, 

.234324 

27. 

,270000 

22. 

.289429 

2 

42, 

,205896 

30, 

,422669 

28. 

390000 

18, 

,639443 

3 

42, 

.38C366 

30, 

,596566 

28. 

,400000 

18, 

,812766 

4 

52, 

.642602 

41  , 

,020917 

51. 

,610000 

29, 

.39923? 

5 

49, 

.639C78 

38, 

,  147348 

4C. 

C40000 

26, 

.655619 

6 

48, 

.794931 

37, 

,317595 

39, 

,  150000 

25, 

,840259 

7 

48, 

,  169940 

36, 

.699251 

41. 

,020000 

25, 

.228561 

8 

48, 

,794931 

37, 

,317595 

38. 

,200000 

25, 

,840259 

9 

49, 

.955463 

38, 

.462943 

34. 

,570000 

26, 

.970424 

10 

48, 

.794S31 

37, 

.317595 

35. 

910000 

25, 

,840259 

11 

48, 

.794931 

37 

.317595 

50. 

,980000 

25, 

.840259 

12 

48, 

,794931 

37, 

.317595 

37. 

,340000 

25, 

.840259 

13 

49, 

,866924 

38, 

.382012 

37. 

320000 

26, 

.895100 

14 

49. 

.789118 

38, 

.300706 

33. 

,690000 

26, 

,812293 

15 

24, 

.544679 

12, 

.757103 

11, 

,320000 

C 

.969526 

16 

22, 

.891935 

11, 

.106287 

11. 

,510000 

-0, 

,679361 

17 

22, 

.891935 

11 

.106287 

12. 

,080000 

-C, 

,679361 

18 

22, 

.891935 

11, 

.106287 

12, 

,490000 

-0, 

.679361 

19 

23, 

.634649 

11, 

.745668 

8. 

,700000 

-0, 

.  143313 

20 

21, 

.821871 

9, 

.921718 

12. 

,930000 

-1, 

.978434 

21 

29. 

.83C509 

18, 

,248583 

18. 

490000 

6, 

,666657 

22 

38, 

.584267 

26 

.911125 

25, 

,  C90000 

15, 

.237882 

23 

34, 

.969159 

23, 

.383477 

23. 

,990000 

11 

.797795 

24 

24, 

.821636 

23, 

.241512 

44. 

7CC000 

11, 

.661383 

25 

25, 

,865046 

24, 

.  165638 

20. 

,250000 

12, 

,466230 

26 

47, 

,221375 

35, 

.535042 

38. 

,490000 

23 

.848708 

27 

43, 

.93C199 

32, 

.312108 

27, 

,270000 

20, 

,694018 

sijfcjfrijtsiijA*****;^^*****  ***  A*:*:*:*:^**:****:**  *«£*  *;***£  ***£:&  :*>!>:#  surest  :£  A;*:*  A  A  #  ri:  »:  A  s  :i 
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TABLE    OF    LS    (INDEPEND)     PREDICTION    INTERVAL    VALUES 
PI     +  PREDICTED  ACTUAL  PI    - 


1 

45, 

.244C64 

34 

.234222 

27. 

269989 

23. 

,224365 

2 

41  , 

.287399 

30 

.422546 

28. 

,389999 

19, 

.557678 

3 

41, 

,473557 

30 

.596436 

28. 

399994 

19, 

.719299 

4 

51, 

,648346 

41 

.020844 

51. 

6C9985 

30, 

.393326 

5 

48, 

,606689 

38 

.147217 

40. 

039993 

27, 

.687729 

6 

47, 

.  776276 

37 

.317474 

39. 

149994 

26, 

.8586  5  8 

7 

47, 

.  151  108 

36 

.699127 

41. 

,019989 

26, 

.247131 

8 

47, 

,776276 

37 

.317474 

38. 

,  199997 

26 

.858658 

9 

48, 

,924454 

38 

.462814 

34. 

569992 

28, 

.001160 

10 

47, 

.776276 

37 

.317474 

35. 

,909988 

26, 

.858658 

11 

47, 

.776276 

37 

.317474 

50, 

,979996 

26, 

.858653 

12 

47 

.776276 

37 

.317474 

37, 

,339996 

26, 

.858658 

13 

48 

.839310 

38 

.381882 

37. 

,319992 

27 

.924438 

14 

48, 

.758224 

38 

.300552 

33. 

689987 

27, 

.842865 

15 

23 

.405258 

12 

.757086 

11, 

,320000 

2, 

.  108913 

16 

21, 

.751877 

11 

.  106257 

11. 

,509999 

0 

.460636 

17 

21, 

.751677 

11 

.  106257 

12. 

C80000 

0, 

.460636 

18 

21 

.751877 

11 

.  106257 

12. 

,490000 

0, 

.460636 

19 

22 

.4333S5 

11 

.745620 

8. 

,700000 

1 

.057642 

20 

20, 

.619278 

9 

.921659 

12. 

,  929999 

-0, 

.775964 

21 

28 

.684082 

16 

.243439 

18, 

,489990 

7 

.812836 

22 

37, 

.507278 

26 

.911041 

25. 

089996 

16 

.314789 

23 

33 

.909637 

23 

.383377 

23, 

,989990 

12 

.857104 

24 

33 

.762192 

23 

.241409 

44. 

,69999  7 

12 

.720620 

25 

24, 

.311674 

24 

.  165512 

20. 

250030 

13 

.519141 

26 

46 

. 174057 

35 

.534973 

38, 

,489990 

24, 

.895374 

27 

42 

.866684 

32 

.312012 

27, 

,269989 

21 

.757324 

*J.  «JU  Jj  -J-   «J.  «J*   -JL   J»  j«  >w   «J*  -J  •    ~£t  «J-   -J*  «J>  «JL»   -J*  «X«  w   v*  •»!■«  «  *  *1*  JL*   «JL>  w  «JU  »i*  JL»  »-*  -J'  *^j   -X-  ■£«  <Jg    J*   *■  '*  -V  *Jg  ~"-*  «X-  >  ■    J—  Of  -  V   -V  'J*  *»-*    -1*  -  "     *^  **-*   **»   ~^*  ^*  ^  -J- 
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MODEL  D 


4  CERs  -  36  OBSERVATIONS 


MULTIPLE  REGRESSION GENRAL 

SELECTION L 

SELECTION  CARD  037  12  L 

LEAST  SQUARES  REGRESSION  COEFICIENTS 


0.470823     0.001158 
STD.  ERROR  OF  ESTIMATE= 


0.804216 


COVARIANCE  MATRIX  OF  THE  ESTIMATES 

0.13244D  00  -0.35830D-04 
-0.3583CD-04   C.  112  14D-07 

ADJUSTED  MULTIPLE  CORRELATION  COEFICIENT  0.772147 


MULTIPLE  REGRESSICN GENRAL 

SELECTION 2 

SELECTION  CARD  038  23341 

LEAST  SQUARES  REGRESSION  COEFICIENTS 


0.005228     0.286414 
STD.  ERROR  OF  ESTIMATE= 


-0.651530 
1.070823 


COVARIANCE  MATRIX  OF  THE  ESTIMATES 

0.26784C  CO  -C.21842D  C2   0.366S9D  02 

-0.21842D  02   C.36286D  C4  -0.21771D  05 

0.36699D  C2  -0.21771D  C5   0.21581D  06 

ADJUSTED  MULTIPLE  CORRELATION  COEFICIENT  0.891135 
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MULTIPLE  REGRESSION GENRAL 

SELECTION 3 

SELECTION  CARD  039  22128 

LEAST  SQUARES  REGRESSION  COEFICIENTS 


-6.785734     O.0C585C 
STD.  ERROR  OF  ESTIMATE= 


6.664107 
4.824125 


COVARIANCE  MATRIX  OF  THE  ESTIMATES 

0.59974D  CI  -C.21718D-C2   0.21987D  01 

-0. 217180-02   0.10359D-C5  -0.15753D-02 

0.21987D  CI  -0.15752D-C2   0.39246D  01 

ADJUSTED  MULTIPLE  CORRELATICN  COEFICIENT  0.844043 


MULTIPLE  REGRESSION GENPAL 

SELECTION 4 

SELECTION  CARD  04C  21422 

LEAST  SQUARES  REGRESSION  COEFICIENTS 


-1.014043     0.005552 
STD.  ERRCP  OF  ESTIMATE= 


0.002013 
2.233319 


COVARIANCE  MATRIX  OF  THE  ESTIMATES 

0.8557CD  CO  -C.54173D-C3  -0.19542D-05 
-0.54172D-C3  C.43047D-C6  -0.36635D-07 
-0.19542D-05  -C.36635D-C7   0. 683600-07 

ADJUSTED  MULTIPLE  CORRELATION  COEFICIENT  0.824401 
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MATRIX  OF  RESIDUALS  FOR  SELECTIONS 


CASE  NO. 

SELECTION  NUMBER 

1 

2 

3 

4 

1 

0.14339 

-0.00057 

-2.94288 

-5.43512 

2 

0.20339 

C. 00005 

-2. 58288 

C. 59742 

3 

0.48339 

C. 00041 

-2.97288 

-0.33258 

4 

-C. 12263 

-0.00005 

0.31065 

3.61752 

5 

-0.57775 

-0.00155 

0.49625 

1.27806 

6 

-0.57848 

-0.00041 

0.19305 

-1.34539 

7 

-1.41931 

C. 00099 

1.17350 

0.52246 

8 

-0.58848 

0.C0030 

-2.30695 

0.6  3461 

9 

0.26940 

-0.00048 

-4.10020 

-2.24420 

10 

-0. 19848 

C.C0064 

-2.75695 

-1.37539 

U 

-0. 15  648 

-0.00050 

10.26305 

I.  12461 

12 

0.23152 

0.00045 

-3.78695 

0.34461 

13 

0.43655 

-0.00109 

-1.64665 

-0.74307 

14 

-0.36260 

-0.00116 

-3.56135 

-0.63627 

15 

0. 83358 

-C. 00115 

1.19070 

-3.01340 

16 

0.39358 

0.00114 

0.78070 

-0.51681 

17 

0.64358 

-C.C0048 

1.07070 

0.34319 

18 

0. 87358 

0.C0151 

1.11070 

0.25319 

19 

-0.3L749 

-0.00056 

-0.44880 

-2.  12091 

20 

0. 12251 

0.00116 

1.12120 

2.08657 

21 

-0.23298 

-0.00123 

0.72584 

-0.66180 

22 

-0.42351 

-0.00064 

-3.30648 

0.47806 

23 

-0.68C65 

-C. 00066 

-2.32003 

3.01800 

24 

-0.87518 

-C. 00025 

15. L3867 

5. 20251 

25 

-0.53170 

-0.00037 

-3.49378 

-0.29083 

26 

-0.48941 

-C. 00036 

2.43806 

-1.28208 

27 

-0.47153 

-0.00033 

-4.35519 

-0.08743 

28 

1.95236 

C. 00212 

5.04403 

4.84444 

29 

-0.  15714 

C. 00077 

-4.77597 

0.  14444 

30 

0.69  80  1 

0.00134 

10.34475 

1.63506 

21 

0.37146 

-C.C0015 

5.04655 

1 .43833 

32 

1.93801 

0.00042 

-7.40525 

-3.41301 

33 

1.56801 

0.00130 

-5.41525 

-1.06301 

34 

-C. 75720 

-C. 00024 

4.3083C 

1. 10306 

35 

-1.30152 

-0.00069 

-1.21455 

-2.53988 

36 

-0.91731 

0.00029 

-0.86380 

-1.56988 

MATRIX  CF  CORRELATIONS  BETWEEN  SELECTIONS 
SELECTION  NO.      I  2  3 


1 

1.000000 

0.441009 

-0. 106147 

-C. 054891 

2 

0.44LC09 

l.COOOOO 

0.084968 

0.278739 

3 

-C. 1C6147 

C.C84968 

1 .000000 

0.555159 

4 

-0.054891 

0.278739 

0.555159 

1.000000 
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TABLE    OF    PREDICTION    INTERVAL    VALUES 
PI    +  PREDICTED  ACTUAL  PI    - 


1 

46, 

.779931 

33. 

,966494 

27. 

270000 

21, 

,153057 

2 

41, 

.496805 

28, 

.721677 

28, 

,390000 

15. 

,946549 

3 

41  , 

,567747 

28, 

.791345 

28. 

400000 

16, 

,014943 

4 

57, 

.000953 

44, 

.309232 

51. 

6  10000 

31, 

,617510 

5 

52, 

.267304 

39, 

.608704 

40. 

,040000 

26, 

.950104 

6 

51. 

.668665 

39, 

,009967 

39. 

150000 

26, 

,351269 

7 

50, 

.944936 

38, 

.280008 

41. 

,020000 

25, 

,615080 

8 

51, 

.668665 

39, 

,009967 

38, 

,200000 

26, 

.351269 

9 

52, 

.574953 

39, 

.918795 

34. 

57C000 

27, 

.262637 

LO 

51. 

.668665 

39, 

.009967 

35, 

,910000 

26, 

.351269 

11 

51, 

.668665 

39, 

,0C9<367 

50. 

980000 

26, 

,351269 

12 

51, 

.668665 

39, 

.009967 

37, 

,340000 

26, 

,351269 

13 

52, 

.836926 

40, 

.183513 

37. 

,320000 

27, 

,530101 

14 

52, 

.655C75 

40, 

.000924 

33. 

690000 

27, 

.346774 

15 

25, 

.302176 

12, 

.369640 

11. 

,320000 

-0. 

,562896 

16 

22, 

.971987 

10, 

.045793 

11, 

,510000 

-2, 

.880401 

17 

22, 

.971987 

10, 

.045793 

12. 

080000 

-2. 

,880401 

18 

22, 

.971987 

10, 

.045793 

12. 

,490000 

-2, 

,880401 

19 

24, 

.486919 

11, 

.580702 

8, 

,700000 

-1. 

.325514 

20 

21 

.968723 

9, 

.C7C958 

12. 

,930000 

—  i 

,826807 

21 

31  , 

.529251 

18, 

.733967 

18. 

,490000 

5, 

.938682 

22 

42, 

.126864 

29, 

.364184 

25. 

,090000 

16. 

,601505 

23 

36, 

.868823 

24, 

.  138036 

22, 

,990000 

11. 

,407248 

24 

26, 

.853816 

24, 

.125882 

44. 

,700000 

11, 

.397948 

25 

36, 

.929294 

24 

.211309 

20, 

250000 

11, 

.483325 

26 

49, 

.691832 

36, 

.953751 

38. 

,490000 

24, 

.215669 

27 

45 

.27C486 

32 

.558521 

27, 

,270000 

19, 

.846556 

28 

60 

.622622 

47 

.915836 

67, 

,590000 

35, 

.208151 

29 

60 

.631846 

47 

.924124 

46, 

,230000 

25, 

.216402 

30 

67 

.494562 

74 

.356561 

93, 

,900000 

61, 

.218559 

31 

79 

.445994 

66 

.331734 

75. 

,310000 

.217575 

32 

78 

.264738 

65 

.180421 

59, 

.300030 

52 

.096105 

33 

78 

.264738 

65 

.  18C421 

62, 

,900000 

52, 

.096105 

34 

87 

.973501 

74 

.452855 

80, 

,040000 

6C, 

.932210 

35 

77 

.611894 

64 

.175077 

59, 

.500000 

50 

.738260 

36 

77 

.615738 

64 

.176842 

62, 

, 170000 

50, 

.737946 

MODEL    D 
JOINT    GENERALIZED    LEAST    SQUARES 
PREDICTION    INTERVAL    VALUES 
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TAELE  OF  LS  CORRELATED  PREDICTION  INTERVAL  VALUES 
PI  +  PREDICTED     ACTUAL        PI  - 


I 

47, 

.496196 

34, 

.597153 

27. 

,  270000 

21, 

,698110 

2 

41, 

,673262 

28, 

.833934 

28. 

,390000 

15, 

,994607 

3 

41, 

,752755 

28, 

.912179 

28. 

4CC000 

16, 

,070603 

4 

57, 

.334560 

44, 

.605451 

51. 

,610000 

31, 

,876342 

5 

52, 

,027124 

39 

.350252 

40. 

,  C40000 

26, 

.673380 

6 

51, 

,422412 

38, 

.745332 

39. 

,  150000 

26, 

,068252 

7 

50, 

.687839 

38, 

.001443 

41, 

,020000 

25, 

.315047 

8 

51, 

,422412 

38, 

.745332 

38. 

200000 

26, 

.06S252 

9 

52, 

.34CE85 

39 

.667734 

34, 

,570000 

26, 

,994582 

10 

51, 

,422412 

38, 

.745332 

35. 

,910000 

26 

.068252 

11 

51, 

.422412 

38 

.745332 

50, 

,  980000 

26, 

.068252 

12 

51, 

.422412 

38 

.745332 

37, 

.340000 

26, 

,068252 

12 

52, 

,614181 

39 

.944989 

37, 

,320000 

27 

.275797 

14 

52. 

.421483 

39 

.760984 

33, 

,690000 

27, 

.090484 

15 

25, 

.116778 

12 

.164660 

11, 

,320000 

-0, 

.787458 

16 

22, 

,551253 

9 

.607907 

11. 

510000 

-3, 

.335440 

17 

22, 

.551253 

9 

.607907 

12, 

,C80000 

-3 , 

.235440 

18 

22, 

.551253 

9 

.607907 

12, 

,490000 

-3 

.235440 

19 

24, 

.283593 

11 

.355201 

8. 

,700000 

-1, 

.573191 

20 

21, 

.514C75 

8, 

.597391 

12. 

,930000 

-4, 

.319294 

21 

21, 

.594248 

18 

.775287 

18, 

,490000 

.956326 

22 

42, 

,729326 

29 

.902177 

25. 

,  C90000 

17, 

,075028 

23 

26, 

.954224 

24 

.  177540 

23. 

,990000 

11, 

.400856 

24 

26, 

.935584 

24 

.163164 

44, 

,700000 

11, 

.390745 

25 

27, 

.018C18 

24 

.245725 

20. 

,250000 

11, 

.473432 

26 

49, 

.791322 

37 

.002960 

38, 

,490000 

24 

.214598 

27 

44, 

.882587 

32 

.133265 

27. 

,270000 

19, 

.382943 

28 

60 

.699060 

47 

.980272 

67. 

,590000 

•3  C 

.261484 

29 

60, 

.707653 

47 

.988864 

46. 

,230000 

35" 

.270076 

30 

88, 

.38323C 

75 

.  122716 

93, 

, 9CC000 

61, 

.862201 

21 

79, 

.428899 

66 

.233432 

75, 

,310000 

53, 

.037964 

32 

78, 

.2C9221 

65 

.039599 

59. 

.300000 

51, 

.869977 

33 

78 

.209221 

65 

.C39599 

62, 

,900000 

51, 

.869977 

34 

89, 

.737422 

76 

.033214 

80, 

,040000 

62 

.339005 

35 

78, 

.311529 

64 

.757111 

59. 

, 500000 

51, 

.202693 

36 

78 

.317571 

64 

.760783 

62, 

, 170000 

51 

.203995 

sjeac****^*^^;*^  £  *******  ric  A***  jfe*-*^**  ■* 
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TAELE  OF  LS 

(INDEPEND)  PRED 

iCTICN  INTE 

RV£L  VALUES 

PI  + 

PRECICTED 

ACTUAL 

PI  - 

1 

45.842146 

34.597580 

27.269989 

23.252997 

2 

40.C2C  11  1 

26  .835066 

28.389999 

17.650009 

3 

40.093536 

28  .913330 

28.399994 

17.733109 

4 

55.642090 

44.605728 

51.609985 

22.569351 

5 

50.345749 

39.251234 

40.039993 

28.256903 

e 

49.722260 

26.746384 

39.149994 

27.770493 

i 

48.965247 

38.002487 

41.019969 

27.039612 

e 

49.72226C 

38.74t384 

38.199997 

27.770493 

c 

5C. 666242 

39.668839 

34.569992 

26.669413 

1C 

49.72226C 

38.746284 

35.909988 

2  7.770493 

11 

49.72226C 

38  .746384 

50.979996 

27.770493 

12 

49.72226C 

36.746264 

37.339996 

2  7.770493 

12 

50  .929285 

39.946075 

37.219992 

28.952350 

14 

5C.747C25 

39.7fc207C 

33.689937 

2  6.7  77100 

15 

23.  195  ItC 

12.  164678 

11.32CC00 

1.134190 

16 

20.629440 

9.608246 

11.509999 

-1.412955 

17 

20.62944C 

9  .60824c 

12. C800Q0 

-1.412938 

18 

20  .62 944 C 

9  .e08246 

12.490000 

-1.412956 

19 

22.365707 

11.355088 

8.700000 

0  .344466 

2C 

19.597290 

8.59761C 

12.929999 

-2.402031 

21 

29.731567 

18.775230 

16.469990 

7.819076 

22 

4C.92527£ 

29.902405 

25.089996 

18.879517 

23 

25.  152029 

24.178421 

23. 969990 

12.204800 

24 

25.122599 

24.164032 

44.699997 

12.195461 

25 

25.2  19786 

24  .246597 

20.250000 

12.273404 

26 

48.CCCC46 

37  .002026 

38.489990 

26.006012 

27 

43.076568 

22.134003 

27.269989 

21  .139423 

2£ 

59.C24C42 

47.9  816  74 

67.589996 

26.929291 

29 

59  .042648 

47.990265 

46.229996 

26.937866 

3C 

66.697797 

75.123077 

93.899994 

63  .348343 

21 

77.8  7  5992 

66  .235046 

75.209993 

54.594086 

22 

76.669571 

65  .041133 

59.299938 

53.412781 

•a  - 

76.669  571 

65.041163 

62.899994 

53.412781 

34 

68.23481£ 

76.036816 

80.029993 

6  3.642  804 

35 

76.781921 

64.759079 

59.500000 

52.736221 

26 

76.7857C7 

64.762741 

62.  169996 

52.739700 
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MODEL  E  -  4  CERs  -  27  OBSERVATIONS 


MULTIPLE  PEGRESSICN ESCCRT 

SELECTION 1 

SELECTION  C/lRC  037  121 

LEAST  SCLARES  PEGRESSICN  CCEFICIEi\TS 

1.2C717C     C. 000825 
STC.  ERRCP  CF  ESTIMATE=        0.49512C 

CCVAPIANCE  MATRIX  CF  THE  ESTIMATES 

0.119660  00  -C.  427920-C4 
-C. 427920-04   C.1656CC-C7 

ADJUSTED  NLLTIPLE  CCRRELATICN  COEFICIENT  0.6C6372 


NULTIPLE  PEGRESSICN ESCCRT 

SELECTION 2 

SELECTION  C£RD  026  122 

LEAST  SCU4RES  REGRESSION  COEFICIENTS 

0.662174    -C. 755852 
STC.  ERRCP  CF  ESTIMATE*        0.912C00 

CCV£PI/iNCE    MATRIX    CF    THE    ESTIMATES 

0.22E57C  CC  -C.2257CD  CI 
-0.2257CC  01   C. 257570  C2 

ADJLSTED    MLLTIPLE    CCRRELATICN    COEFICIENT    0.894681 


99 


NULTIPLE  REGRESSICN ESCCRT 

SELECTICN 3 

SELECTICN  C/^RC  C39  ^212£ 

LEAST  SCIAPES  REGRESSICN  CCEFICIENTS 


-4.1883C1     C. 004682 
STC.  EPRCR  QF  ESTIMATE= 


7.148997 
4.414830 


CCVARIANCE  MATRIX  CF  THE  ESTIMATES 

0.13266C  C2  -C.57465C-C2   0.47872D  01 

-0.574650-C2   C.27811C-C5  -0.299CcD-02 

0.47872D  CI  -C.299C6D-C2   C.61C73D  01 

ADJUSTED  NLLTIPLE  CCRRELATICN  COEFICIENT  0.687568 


MJLTIPLE  REGRESSICN ESCCRT 

SELECTICN 4 

SELECTICN  CAPC  C4C  21423 

LEAS1  SCLAPES  REGRESSION  CCEFICIENTS 


-2.627433  •   0.0C7357 
STC.  ERRCR  OF  ESTIMATE= 


0.0C1560 
2.C51149 


CCVARIANCE  MATRIX  CF  THE  ESTIMATES 

0.26712C  01  -C.23852C-C2  -C.S75C3D-04 
-C.23E52D-C2  C.22C61C-C5  0.149670-07 
-0.97SC2C-C4   C.14967D-C7   0. 139220-06 

ACJLSTED  VLLTIPLE  CCRRELATICN  COEFICIENT  0.593167 
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MATRIX  CF  RESIDLALS  FCP  SELECTICNS 


CASE    NG. 

SELECTION    NUMBER 

1 

2 

3 

4 

i 

0.26221 

-0.01455 

-2.19580 

-4.25904 

2 

0,42221 

C.C3284 

-1.83580 

0.47754 

3 

0.7G221 

C.C3799 

-2.225SC 

-0.45246 

4 

0.  19227 

C  .03229 

0.91279 

4.47934 

5 

-C.24C5C 

-C.C6G26 

1.17315 

0.671^0 

6 

-0.24290 

C.C1815 

C.3606C 

-1.86567 

7 

-1.11C41 

C. 10234 

1.75111 

-0.05373 

8 

-C.2539C 

CO  7  53  5 

-1.63940 

0.  11433 

9 

0.fcL90C 

0.C2679 

-3.33008 

-2.8258C 

10 

0.13  61C 

0.10293 

-2.08940 

-1.89567 

11 

C.  176  1C 

C .01108 

1C.9306C 

0.60433 

12 

0.56610 

0 .08756 

-3.  U94C 

-0.  17567 

13 

C. 79849 

-C.C1591 

-0.33331 

-1.33730 

14 

-C.0C667 

-0.02876 

-2.31903 

-1.25765 

15 

0.53C94 

-C. 19377 

0.11191 

-1  .85765 

16 

C.C9C94 

-C.1C238 

-0.29809 

0.05043 

17 

C.34C94 

-C.  16972 

-0.00309 

0.91048 

18 

C.57C94 

-C. 03759 

0.03191 

0.82043 

19 

-C. 59677 

C.CC624 

-1.44582 

-0.97135 

20 

-0.  15677 

C .  12932 

0.12418 

2.61599 

21 

-C. 35512 

-C.CC047 

0.279C4 

0.34926 

22 

-0.26608 

-0. 02409 

-2.77434 

1. 10732 

23 

-0.53557 

-0.02506 

-2.33112 

2.45983 

24 

-0.73743 

C.CC409 

15.o0138 

4.59384 

25 

-C. 39496 

C.CC603 

-3.03  40  7 

-0.89409 

26 

-0.22164 

C  .00365 

2.52152 

-0.66404 

27 

-C. 26874 

C.C2137 

-4.21916 

-0.74479 

NATRIX  CF  CORRELATIONS  EETV\EEN  SELECTICNS 
SELECTICN  NC.      1  2  3 


l.OOCCOC 
-0.255531 
-C. 280213 
-C. 317896 


0.255531 

l.CCCCOO 

•0.043023 

C.C7C142 


-0.280213 

-C. 043023 

l.CCOGOC 

C. 490192 


■0.317896 
C. 070142 
C.49C192 

1.000000 
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TABLE  CF  PREDICTICN  INTERVAL  VALUES 
FI  +  PREDICTED     ACTUAL        PI  - 


1 

42.744578 

31.861400 

27.270000 

2C. 978222 

2 

38.1443C2 

27.352d79 

28.390000 

16  .  56143b 

3 

28 . 566  898 

27.773446 

28.400000 

16.979995 

X 

53.43C6<5 

42.696632 

51.610000 

3  1.962670 

5 

43.9C7651 

33.248580 

40. 040000 

27.589510 

6 

48.522766 

37.87742C 

39.  150000 

27.232095 

7 

48.297413 

37.651677 

41.020000 

27.005941 

8 

48.522766 

37.87743C 

28.2GC000 

27,232095 

9 

49.023746 

38.363142 

34.57C000 

27.702538 

1C 

48.522766 

37.87743C 

35.91C000 

27.232095 

LI 

48. 522766 

37.877430 

5C.98C000 

27.232095 

12 

48.522766 

37  .877430 

27.340000 

27.232095 

13 

4  9.2  89304 

38.62C256 

37.320000 

27.951207 

14 

4  9  . 2  9  C  2  L  6 

33  .62CL12 

33.69C000 

27.949909 

15 

24.6380ol 

13.611625 

11.320000 

2.585188 

16 

22.729269 

LI. 706695 

11. 510000 

C.  684520 

17 

22.729269 

LI. 706895 

12.C80000 

0.684520 

18 

22.7292e9 

11  .706895 

12.490000 

0.684520 

19 

21.643423 

10.869270 

8.7CC000 

-0.  104882 

2C 

19.679577 

S. 711857 

12.930000 

-2  .255862 

21 

27.8534C1 

L7. 156160 

18.490000 

6.453919 

22 

38.739543 

27.946L16 

25.C9C000 

17.  L52688 

23 

24.119218 

2^.422719 

23.990000 

12.726120 

24 

34. 1 7S  6  IS 

23.468163 

44.7CC000 

12.796707 

25 

34.  14C610 

23.450805 

2C. 250000 

12.760999 

26 

46.75  2  112 

3O.022949 

38.49C000 

25.293766 

27 

42.864222 

32.193560 

27.270000 

21.522899 

MODEL  E 
JOINT  GERERALIZED  LEAST  SQUARES 
PREDICTION  INTERVAL  VALUES 
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TAELE    OF    LS    CORRELATED    PREDICTION    INTERVAL    VALUES 
PI     +  PREDICTED  ACTUAL  PI    - 


1 

42, 

,996416 

31, 

.976923 

27. 

270000 

20, 

,957429 

2 

38, 

.402254 

27 

.510643 

28. 

,390000 

16, 

,618931 

3 

38, 

.815021 

27 

.921L90 

28. 

,400000 

17, 

,027360 

4 

,056183 

42 

.264582 

51. 

610000 

31, 

.472982 

5 

48, 

,544981 

37, 

.866001 

40, 

,040000 

27. 

,  187021 

6 

48, 

.133818 

37 

.468809 

39, 

,150000 

26, 

.803801 

7 

47, 

,942283 

37 

.272837 

41. 

,  020000 

26. 

,603392 

8 

48, 

.133818 

37, 

.468809 

38, 

,200000 

26, 

.803301 

9 

48, 

,645678 

37 

.966762 

34. 

,570000 

27, 

,287847 

10 

48, 

.  132818 

37 

.468809 

35, 

,910000 

26, 

.803801 

11 

48, 

.133818 

37 

.468809 

50, 

,980000 

26, 

.803801 

12 

48, 

,  133818 

37 

.468809 

37. 

,340000 

26, 

,803801 

13 

48, 

.908250 

38 

.222020 

37. 

,320000 

27, 

,535791 

14 

48, 

,921407 

38 

.233249 

33, 

,690000 

27, 

.545092 

15 

24, 

.942192 

13 

.883649 

11, 

,320000 

2, 

,825107 

16 

23, 

.050856 

11 

.994706 

11, 

,510000 

0, 

.938556 

17 

23, 

.05C856 

11 

.994706 

12, 

,080000 

0, 

.938556 

18 

23, 

.050856 

11 

.994706 

12, 

,490000 

0, 

.938556 

19 

22, 

.225079 

11 

.221966 

8, 

,700000 

0 

.218853 

2C 

20, 

.082272 

9 

.084766 

12. 

,  930000 

-1, 

.912739 

21 

28, 

.22C243 

17 

.501517 

18. 

,490000 

6, 

,782790 

22 

39, 

.242649 

28 

.342673 

25, 

,090000 

17 

.442697 

23 

24, 

.622194 

23 

.856676 

23, 

,990000 

13, 

.091159 

24 

24, 

.699258 

23 

.941518 

44, 

,700000 

13, 

.  183779 

25 

24, 

,658616 

23 

.9  0292  6 

20, 

,250000 

13 

.147237 

26 

46, 

.575587 

35 

.767393 

38, 

,490000 

24, 

.959200 

27 

42, 

.711111 

31 

.986247 

27, 

.270000 

21 

,2ol384 

-  ■-       -1-*    —   -        <-     -.' -       -  •    Uj    •    «JU    U*    ~—r    -J  *    aX.     *  »    kjg    -J--       ~*     tf  -    V*     -X-     J I  '    - — •    ••  *     J* 

«i «    >!#■-'*   -J*    — ^    V-»  -J^  %*#  **£    *i*  *JL>  O'  >Ju   o<   C-  »:»  w     .'*  vL*     -4-   •.'-*    ~*j    "'-    v~"  »**  »J-»    -  ■* 


rt  ri:  J:  jis  -i  is  •&  A  -it  xt  Jis  -^  •?*  i:  A  5fe  *  ~  '^  "^  "  -:  5't  ^ 
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TAeLE  CF  LS 

(INDEPEND)  PREDICTION  INTERVAL  VALUES 

PI 

+ 

PREDICTED 

ACTUAL 

PI  - 

1 

42, 

,668549 

31.974792 

27.269939 

21.281021 

2 

38, 

.081787 

27.5CS209 

28.289999 

16.924616 

-a 

28, 

,5  1799 C 

27.918762 

28.399994 

17.319519 

4 

52. 

,661728 

42.262314 

51.6C9985 

31.662885 

5 

48, 

.1C5C57 

37.663342 

40.C29992 

27.621cl3 

6 

47, 

.7CCC72 

37.466171 

29.149994 

27.232254 

7 

47, 

.5C5875 

37.^70248 

41.C19939 

27.C34607 

8 

47, 

.7CCC73 

37  .466171 

38.199997 

27.232254 

9 

48, 

.2C72C6 

37.964081 

34.569992 

27.720340 

LC 

47, 

.7CCC72 

37  .466171 

35.9C9933 

2  7.22  22  54 

11 

47, 

.7CCC72 

37  .46fcl71 

50. 9799^6 

27.222254 

12 

47, 

.7CCC72 

37.466171 

27.329996 

27.222254 

12 

46, 

.466902 

38.219299 

27.319992 

27.969681 

14 

48, 

,479965 

38.23C545 

33.639987 

27.981110 

15 

24, 

,218436 

13.882758 

11.320000 

3.447C76 

16 

22, 

.433441 

11  .993677 

11.5C9999 

1.553909 

17 

22 , 

,43j441 

11  ,99^c77 

12.C30000 

1.553909 

16 

22  , 

.42344  1 

11.993677 

12.490000 

1.553909 

19 

21  , 

,  51 6C36 

11.221060 

8.700000 

C. 924072 

2C 

19, 

,37  / 1  82 

9.C83716 

12.929999 

-1.20976i 

21 

27, 

,562980 

17.500229 

18.489990 

7.436467 

22 

28 , 

.695C3S 

28.34C729 

25.08CJ996 

17  .986404 

23 

24, 

,C7£553 

23.654492 

23.989990 

12.632421 

24 

24, 

.152678 

23.939346 

44.699997 

13.725811 

25 

24 , 

,  112C3C 

23.900742 

20.250000 

13.689447 

26 

46, 

.C4  6  14  2 

35.765503 

26.489990 

25.464848 

27 

42. 

.136733 

31  .964024 

27.269989 

21.829300 
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APPENDIX  D 
COMPUTER  PROGRAMS 

A.   JOINT  GENERALIZED  LEAST  SQUARES  PROGRAM 
1 .   Program  Description 

The  joint  generalized  least  squares  computer  program 
is  designed  to  read  in  a  set  of  data  and  then  carry  out  joint 
generalized  least  squares  computations  as  specified  on  a  con- 
trol card,  and  selection  cards.   The  program  uses  a  number 
of  IBM  360  Scientific  Subroutine  packages  to  do  the  calcula- 
tions.  It  is  programmed  using  FORTRAN  IV  language  much  like 
the  Multiple  Regression  program  in  the  scientific  subroutine 
package  (REGRE).   The  program  conducts  a  number  of  standard 
multiple  linear  regressions  (one  for  each  subsystem  CER) . 
The  residuals  from  the  single  regresions  are  saved  and  uti- 
lized to  compute  a  matrix  of  correlations  between  selections. 
Both  of  these  items  are  listed  in  the  output.   Next  the 
program  does  the  joint  generalized  least  squares  calcula- 
tions; initially  S,  the  matrix  of  mean  squares  and  products 
of  the  least  squares  residuals,  is  computed,  next  b.,  the 
joint  parameter  estimate  vector  is  computed  along  with  V(b.), 
the  covariance  matrix  of  the  joint  estimates.   Finally,  a 
set  of  joint  generalized  least  squares  prediction  interval 
values  is  calculated.   The  program  was  designed  to  be  easy 
to  use  for  analysis  utilizing  the  joint  generalized  least 
squares  method.   Double  precision  was  used  to  improve  accuracy 
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by  reducing  roundoff  errors.   The  program  will  handle  up 
to  thirty-six  observations  and  nine  subsystem  relations. 
Each  subsystem  CER  may  have  up  to  three  explanatory  variables 
The  program  takes  280  K  of  core  storage  space  to  run  with 
a  maximum  run  time  of  about  thirty  seconds.   It  should  be 
noted  that  larger  dimensional  problems  may  be  possible  but 
computer  capacity  rapidly  becomes  a  constraint  to  be  con- 
cerned with.   The  scientific  subroutine  ARRAY  was  used  to 
allow  variable  dimensioning  so  that  those  subroutines  that 
deal  with  matrices  can  operate  on  any  size  array,  limited 
only  by  the  maximum  program  dimensions . 
2 .   JGLS  Program  Users  Guide 

The  program  listed  herein  is  straightforward  to  use. 
The  user  must  supply  as  input  a  data  deck  of  independent  and 
dependent  variables.   The  data  array  is  presently  programed 
to  store  up  to  forty-one  variables  and  up  to  thirty-six 
observations  for  each  variable.   The  user  should  supply  the 
appropriate  read  statements  in  the  program  to  accommodate 
the  data  format  he  is  using.   After  the  data  is  read  in  any 
transformations  desired  may  be  made. 

a.   Control  Card 

One  control  card  is  required  for  each  program 
and  is  read  before  the  data  deck.   This  card  is  prepared 
as  follows: 
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COLUMNS 


CONTENTS 


SAMPLE 


1-6 
7-11 
12-13 
14-15 


Problem  Name  Escort 

Number  of  Observations  36 

Number  of  Variables  4l 

Number  of  Selection  Cards  04 


b.   Selection  Card 

The  selection  card  is  used  to  specify  the  type 
of  transformation  desired,  the  dependent  variable  and  a  set 
of  independent  variables  in  a  subsystem  multiple  linear 
regression  analysis.   One  selection  card  is  required  for 
each  subsystem  in  the  model.   The  variables  are  designated 
as  columns  in  the  data  array.   The  selection  card  is  prepared 
as  follows: 

COLUMNS         CONTENTS 


1-2 

3-4 

5-6 

7-8 
9-10 
11-12 


Transformation  Code 

00  Linear 

01  Log  Linear 

Dependent  Variable 
Designated  for  the 
Regression 

Number  of  Independent 
Variables   Included  in 
Regression 

st 

1  Independent  Variable 

Included 

2  Independent  Variable 
Included 

rd 

3  Independent  Variable 

Included 


SAMPLE 


00 


37 

01 
21 
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The  input  format  3612  is  used  for  the  selection  card. 
c.   Deck  Setup 

The  deck  setup  to  run  the  program  for  a  number 
of  selections  is  as  follows: 

1.  Main  Program 

2.  Control  Card 

3.  Data  deck  (Read  statements  must  be  supplied  in  the 
main  program. 

4.  Selection  Cards  (One  selection  card  for  each  subsystem 
multiple  linear  regression) . 

The  joint  generalized  least  squares  program  is  listed  in  the 

following  section  of  this  Appendix. 
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